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Abstract 

Analytical expressions for the saturation density of asymmetric nuclear matter as well as its 
binding energy and incompressibility at saturation density are given up to 4th-order in the isospin 
asymmetry 5 = {pn — Pp)/ p using 11 characteristic parameters defined by the density derivatives 
of the binding energy per nucleon of symmetric nuclear matter, the symmetry energy Esjra{p) 
and the 4th-order symmetry energy £'sym,4(/') at the normal nuclear density po- Using an isospin- 
and momentum-dependent modified Gogny (MDI) interaction and the Skyrme-Hartree-Fock (SHF) 
approach with 63 popular Skyrme interactions, we have systematically studied the isospin depen- 
dence of the saturation properties of asymmetric nuclear matter, particularly the incompressibility 
Kse,t{5) = + i^sat,2^^ + -f^sat,4'5^ + 0{5^) at Saturation density. Our results show that the 
magnitude of the higher-order Xsat,4 parameter is generally small compared to that of the i^^sat,2 
parameter. The latter essentially characterizes the isospin dependence of the incompressibility at 
saturation density and can be expressed as i^sat,2 = -f^sym — 6L — where L and Ksym rep- 

resent, respectively, the slope and curvature parameters of the symmetry energy at po while Jq 
is the third-order derivative parameter of symmetric nuclear matter at po- Furthermore, we have 
constructed a phenomenological modified Skyrme-like (MSL) model which can reasonably describe 
the general properties of symmetric nuclear matter and the symmetry energy predicted by both 
the MDI model and the SHF approach. The results indicate that the higher-order Jq contribution 
to i<rsat,2 generally cannot be neglected. In addition, it is found that there exists a nicely linear cor- 
relation between iTsym and L as well as between Jq/Kq and Kq. These correlations together with 
the empirical constraints on Kq, L, -Esym(po) and the nucleon effective mass lead to an estimate of 
i^sat,2 = -370 ± 120 MeV. 

PACS numbers: 21.65.Mn, 21.65.Ef, 21.30.Fe 
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I. INTRODUCTION 



The equation of state (EOS) of nuclear matter is one of central issues in nuclear physics. 
For a cold nuclear matter, the EOS is usually defined as the binding energy per nucleon as a 
function of the density from which information on other thermodynamic properties of nuclear 
matter, such as its pressure and incompressibility can be obtained. With the establishment 
or construction of many radioactive beam facilities around the world, such as the Cooling 



Storage Ring (CSR) 
at RIKEN in Japan 



factory 
J], and 



acility at HIRFL in China the Radioactive Ion Beam (RIB) '. 
2|, the FAIR/GSI in Germany U, SPIRAL2/GANIL in France 
the Facility for Rare Isotope Beams (FRIB) in USA |5|, it is possible in terrestrial laboratories 
to explore the EOS of an isospin asymmetric nuclear matter under the extreme condition of 
large isospin asymmetry, especially the density dependence of the nuclear symmetry energy. 
This knowledge, especially the latter, is important for understanding not only the structure 
of radioactive nuclei, the reaction dynamics induced by rare isotopes, and the liquid-gas 
phase transition in asymmetric nuclear matter, but also many critical issues in astrophysics 
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For symmetric nuclear matter with equal fractions of neutrons and protons, its EOS is 
relatively well-determined after about more than 30 years of studies by the nuclear physics 
community. In particular, the incompressibility of symmetric nuclear matter at its saturation 



density oq has been determined to be 240 

n n n 

monopole resonances (GMR) pj, 1171, 118|: 
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20 MeV 
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Tom analyses of the nuclear giant 
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241 ] ■ and its EOS at densities 



of 2po < P < 5po has also been constrained by measurements of collective flows in nucleus 
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261 ] in relativistic nucleus- 



nucleus collisions js] and of subthreshold kaon production 
nucleus collisions. On the other hand, the EOS of asymmetric nuclear matter, especially 
the density dependence of the nuclear symmetry energy, is largely unknown. Although 
the nuclear symmetry energ y a t po is known to be around 30 MeV from the empirical 
liquid-drop mass formula 27l.j28l. its values at other densities, especially at supra-saturation 
densities, are poorly known p, LZl]. Various microscopic and phenomenological rno dels , such 



as the relativistic Dirac-Brueckner-Hartree-Fock (DBHF 
the non-relativistic Brueckner-Hartree-Fock (BHF) 



)BHF) [29|, l30 
3, 37, y, 3^ 
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35| and 



approach, the relativistic 



mean-field (RMF) model based on nucleon-meson interactions [l^, 
relativistic mean-field model based on Skyrme-like interactions 
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44 



45|, 



42], and the non 
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49, 
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5l| . have been used to study the isospin-dependent properties of asymmetric nuclear 



matter, such as the nuclear symmetry energy, the nuclear symmetry potential, the isospin- 
splitting of the nucleon effective masses, etc., but the predicted results vary widely. In fact, 



even the sign of the symmetry energy above 3po is still uncertain [52|, |53|]. The theoretical 
uncertainties are mainly due to the lack of knowledge about the isospin dependence of in- 
medium nuclear effective interactions and the limitations in the techniques for solving the 
nuclear many-body problem. 

Heavy-ion collisions, especially those induced by neutron-rich nuclei, provide a unique tool 
to investigate the EOS of asymmetric nuclear matter, especially the density dependence of 
the nuclear symmetry energy. During the last decade, significant progress has indeed been 



made both experimentally and theoretically on constraining the 
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be 
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havior o 
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57 



th e sy mmetry 
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59, 



energy at subsaturation density using heavy-ion reactions 
Ref. 15|] for the most recent review). More recently, the IBUU04 transport mode^ 



60| (See 



analysis 



of the FOPI data on the Ti~ jn^ ratio in central heavy-ion collisions at SIS/GSI 6l| energies 
suggests a very soft symmetry energy at the suprasaturation densities 621]. Information 
on the density dependence of the nuclear symmetry energy has also been obtained from 
the structure of finite nuclei and their excitations, such as the mass data 63 L neutron 



skin in heavy nuclei 6J], giant dipole resonances 65|], pygmy dipole resonance 66|, and so 
on. These studies have significantly improved our understanding of the EOS of asymmetric 
nuclear matter. 

The incompressibility of asymmetric nuclear matter at its saturation density is a basic 
quantity to characterize its EOS. Since this quantity is largely undetermined, any constraint 
imposed on it would be extremely important. In the present work, we study the isospin 
dependence of the properties of asymmetric nuclear matter, including the saturation density 
as well as the binding energy and incompressibility at saturation density. In particular, we 
derive analytical expressions for these quantities up to the 4th-order in the isospin asymmetry 
<^ = (Pn — Pp) / P and investigate the higher-order isospin asymmetry effects on the properties 
of asymmetric nuclear matter. For the incompressibility of an asymmetric nuclear matter 
at its saturation density, it can be written as K^.^x.{S) = -ft'o + -^'sat,2^^ + -^sat,4^'^ + 0(5®) with 
-^sat,2 = -^'^sym — 6L — -^L, whcrc L and /'^sym represent, respectively, the slope and curvature 
parameters of the symmetry energy at po while Jq is the 3rd-order derivative parameter of 
symmetric nuclear matter at po. Therefore, the higher-order effects on -ft'sat(5) also include 
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the 3rd-order density derivative parameter Jq. Our results indicate that higher-order isospin 
asymmetry effects on the incompressibihty are usually not important but the higher-order Jq 
contribution generally cannot be neglected. Furthermore, we construct a phenomenological 
modified Skyrme-like (MSL) model which can reasonably describe the general properties 
of symmetric nuclear matter and the symmetry energy predicted by both the MDI model 
and the Skyrme-Hartree-Fock (SHF) approach. We find that there exists a nicely linear 
correlation between -K'sym and L as well as between Jq/Kq and Kq. These correlations 
together with the empirical constraints on i^'o, L and -EsymlPo) lead to an estimate of i^^sat,2 = 
-370 ± 120 MeV. 

The paper is organized as follows. In Section [ITl we discuss the general properties of 
asymmetric nuclear matter, and then give analytical expressions up to 4th-order terms in 
isospin asymmetry 6 for the saturation density of asymmetric nuclear matter as well as 
its binding energy and incompressibihty at saturation density. We then briefly introduce in 
Section UTTl the three models used in the present paper, i.e., the isospin- and momentum- 
dependent MDI model, the SHF approach, and the phenomenological MSL model. The 
results and discussions are presented in Section llVi A summary is then given in Section 
IVl For completeness, derivations of some important formula shown in Section [III are briefly 
described in the Appendix. 

II. SATURATION PROPERTIES OF ASYMMETRIC NUCLEAR MATTER 
A. Equation of state of asymmetric nuclear matter 

The EOS of isospin asymmetric nuclear matter, given by its binding energy per nucleon, 
can be generally expressed as a power series in the isospin asymmetry 6 = {pn — Pp) / P, where 
P = Pn + Pp is, the baryon density with p„, and pp denoting the neutron and proton densities, 
respectively. To the 4th-order in isospin asymmetry, it is written as 

E{p, 6) = Eo{p) + ^sym(p)5' + ^sym,4(p)5' + 0(5^), (l) 
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where Eo{p) = E{p,S = 0) is the binding energy per nucleon of symmetric nuclear matter, 
and 

Esym[p) - ^ 1 5=0 [^) 

^sym,4lPj = ^ -g^^ |5=0- (-3) 

In the above, Esymip) is the so-called nuclear symmetry energy and -Esym,4(p) is the 4th- 
order coefficient, which is called here the 4th-order nuclear symmetry energy. The absence 
of odd-order terms in 6 in Eq. ([T]) is due to the exchange symmetry between protons and 
neutrons in nuclear matter when one neglects the Coulomb interaction and assumes the 
charge symmetry of nuclear forces. The nuclear symmetry energy -E'sym(p) thus corresponds 
to the lowest-order coefficient. The higher-order (including 4th-order) coefficients in 6 are 
usually very small and neglected, e.g., the magnitude of the 5^ term at the normal nuclear 
density po (the saturation density of symmetric nuclear matter) is estimated to be less than 



1 MeV in microscopic many-body approaches [67|, 168|, 169| and also in phenomenological 
models as shown later. Neglecting the contribution from higher-order terms in Eq. ([T]) leads 
to the well-known empirical parabolic law for the EOS of asymmetric nuclear matter, which 
has been verified by all many-body theories to date, at least for densities up to moderate 



values 



15|. As a good approximation, the density-dependent symmetry energy i?sym(p) can 



thus be extracted from the parabolic approximation of -E'sym(p) ~ E{p,6 = 1) — E{p,6 = 0). 
It should be mentioned that the possible presence of the higher-order terms in 6 at supra- 
normal densities can significantly modify the proton fraction in /^-equilibrium neutron-star 
matter and the critical density for the direct Urea process which can lead to faster cooling 



of neutron stars 



70 



71| . In addition, a recent study indicates that the higher-order 



terms in 6 are very important for determining the transition density and pressure at the 
inner edge separating the liquid core from the solid crust of neutron stars where the matter 
is extremely neutron-rich. 

Around the nuclear matter saturation density po, the binding energy per nucleon in 
symmetric nuclear matter Eq{p) can be expanded, e.g., up to 4th-order in density, as 

Eo{p) = Eoipo) + LoX+^x' 

+ §X^ + §X^ + 0(x^), (4) 
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where x is a dimensionless variable characterizing the deviations of the density from the 
saturation density po of symmetric nuclear matter and it is conventionally defined as 



The first term Eo{po) on the right-hand-side (r.h.s) of Eq. (jll) is the binding energy per 
nucleon in symmetric nuclear matter at the saturation density po and the coefficients of 
other terms are 

_ „ dEo{p) 
^ (fEo{p) 



Jo = 27po — — |p=po, (8) 

. d'^Enip) , , , 

^0 = 81p^^^|,=,o. (9) 

Obviously, we have Lq = according to the definition of the saturation density po of sym- 
metric nuclear matter and thus the second term on the r.h.s of Eq. (j4]) should vanish. 
The coefficient Kq is the so-called incompressibility coefficient of symmetric nuclear matter 
and it characterizes the curvature of Eo{p) at po. The coefficients Jo and Iq correspond to 
higher-order contributions and here we call them as 3rd-order and 4th-order incompress- 
ibility coefficients of symmetric nuclear matter, respectively. In the literature, one usually 
neglects the higher-order terms in Eq. (jlj) around the saturation density po and obtain the 
following parabolic approximation to the EOS of symmetric nuclear matter: 

E,ip) = Eo{po) + ^x' + Oix'). (10) 

Around the normal nuclear density po, the nuclear symmetry energy i?sym(p) can be 
similarly expanded, e.g., up to 4th-order in x, as 

Kg m 

^sym(p) = ^sym(po) + + ^pX^ 

+ + (11) 

where L, -K'sym, Jsym and Jgym are the slope parameter, curvature parameter, 3rd-order 
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coefficient, and 4th-order coefficient of the nuclear symmetry energy at po, i.e. 

'sym 

dp 



3P„%^|,.„, (12) 



J^sym — ^Po o 9 \p=Poy \^^) 



J _ 27 3 ^^-^sym(p) | 
-'sym — ^'Pq lp=po' \^^) 



J _ Q1 4 '^^-^sym(p) I /.f-x 

Jsym — OJ-Po ^^^^i lp=P0- V-'-'Jj 

The coefficients L, /isym, Jsym and /gym characterize the density dependence of the nuclear 
symmetry energy around the normal nuclear density po? and thus carry important informa- 
tion on the properties of nuclear symmetry energy at both high and low densities. 

We can also expand the 4th-order nuclear symmetry energy -E'sym,4(p) around the normal 
nuclear density po up to 4th-order in x as 

-E'sym,4(p) = -£'sym,4(Po) + L 

sym, 4 

+ + (16) 

where L^yra^^ -^sym,4, <^sym,4 and /sym,4 are the slope parameter, curvature parameter, 
3rd-order coefficient, and 4th-order coefficient of the 4th-order nuclear symmetry energy 
-Esym,4(p) at Po, i.e., 

J _ ^ dE,y^^i{p) 

-t^sym,4 — JPO 1 lp=Po' \^') 



dp 

^2 '^^-£'sym,4(p) 



— 0^2 " -^sym,4V^; , /^„x 
-fi-sym,4 — »Po j o IP=PO' \^°) 



<>'sym,4 — ^'Po ^^3 lp=PO) K^^J 

sym. 



^ C?'^£'sylYl,4 (p) I ^ K 

-'sym,4 = 81po -J-J \p=po- (20) 



In above Taylor expansions, we have kept all terms up to 4th-order in 6 or x- The 14 

coefficients, namely, -Eo(Po), Kq, Jq, Jq, -Esym(Po), L, Ksym, ^sym, ^sym, -Esym,4(Po), -^^sym,4, 

-f^sym,4, «/sym,4 and /sym,4, are wcll-defined, and they characterize the EOS of an asymmetric 
nuclear matter and its density dependence at the normal nuclear density po. Among these 
parameters, Eq{pq), Kq, Esym{Po), L, and Ksym have been extensively studied in the literature 
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and significant progress has been made over past few decades. Based on tliese well-defined 
quantities, we investigate in the present paper to what extend they can provide reliable 
information on the isospin dependence of the saturation properties of asymmetric nuclear 
matter as well as its properties at both high and low densities. 



B. Saturation density 

One of the basic quantities for describing an asjmimetric nuclear matter is its saturation 
density Psat(^) which is generally a function of the isospin asymmetry 6. According to the 
definition of the saturation density Psat(5) of asymmetric nuclear matter, i.e., 

dE{p,6). 



dp 



IP=Psat 



(5) = 0, 



(21) 



it can be shown that up to 4th-order in 6, the saturation density can be expressed as (see 
Appendix ^ 



Psat(5) 



'■0 



6^ + 0{6' 



6\ 



Po- 



Therefore, Psat(^) can be written as 



(22) 



Psat(5) = Po + Psat,2^^ + Psat,45'^ + 0(5^) 



(23) 



with 



Psat,2 
Psat,4 



3L 

-To'' 



-"■0 



3-Z^sym,4 SJqL^ 







2Kl 



Po 



(24) 
(25) 



which reflects the shift of the saturation density of asymmetric nuclear matter relative to 
that of symmetric nuclear matter due to the finite isospin asymmetry. 



C. Binding energy at saturation density 

Another basic quantity of asymmetric nuclear matter is the binding energy per nucleon 
at saturation density, i.e., Esa_t{S), and it is generally a function of the isospin asymmetry S. 

9 



According to Eq. ([T]) and Eq. (l23i) . Esa.t{S) can be expressed up to 4th-order in S as 



Esat(5) = ^O(PO) + 
+ 



Kq ( Psat - PO 



2 



-t^symlPoj + 



2 X 9 V Po 

2 , 771 \ r4 



-J V Po /. 
z X y Aq 

-T7i^5' + i^sy.,4(po)5' + 0(5^) 
6 Kq 

= ^o(Po) + ^sym(Po)^^ + ^£'sym,4(Po) " 

+0(5^), (26) 

Therefore, E'satl'^) can be written as 

£^sat(5) = £^o(po) + E,,t,2S^ + E^^^aS^ + 0{S'), (27) 

with 

-E^sat,2 = -£^sym(Po) (28) 

^2 

-E'sat,4 = -E'sym,4(Po) " 777:^ • (29) 

ZKq 

The binding energy of asymmetric nuclear matter at saturation density is thus shifted from 
that of symmetry nuclear matter as a result of the shift of the saturation density when the 
isospin asymmetry is finite. 

D. Incompressibility at saturation density 

The incompressibility of asymmetric nuclear matter is an important quantity to describe 
its properties. It depends on the density and isospin asymmetry of the asymmetric nuclear 
matter and is characterized by the incompressibility coefficient defined as 

dP{p,6) 



Kip, 5) = 9 



dp 

IB^.V^. (30) 
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where P(p, S) is the pressure of asymmetric nuclear matter and it can be expressed as 

(31) 

Conventionally, the incompressibility coefficient is defined at the saturation density where 
P(p, 5) = 0. It is also called the isobaric incompressibility coefficient 7J], and is given by 

i^sat(5)=9pL^^|^lp=P..- (32) 

The isobaric incompressibility coefficient -R'sat(^) thus only depends on the isospin asymmetry 
6. One can show (see Appendix □) that up to 4th-order in 6, the isobaric incompressibility 
coefficient -ft'satl'^) can be expressed as 



I / TV- i^T ^0-^sym,4 9L JgymL 

+ l-"-sym,4 - 0-t^sym,4 77 H -jp 77 

Kq Kq Kq 

+0{6'), (33) 
which can be further written as 

Ks.t{S) = Ko + K,,t,2S^ + K,^,^,5^ + 0(^6) (34) 



with 



-^sat,2 — -^sym — 6L — J^^^ (35) 







<^0-^sym,4 9L^ JsjmL 



-'^sat,4 — -'^sym,4 U-'^sym,4 ' 

Kn An K( 



^0 



The coefficient A'sat,2 and i^sat,4 refiect the isospin dependence of the isobaric incompress- 
ibility of asymmetric nuclear matter. 

If we use the parabolic approximation for the EOS of symmetric nuclear matter, i.e., Eq. 
( fTOl) . then the -R'sat,2 parameter is reduced to 

-^asy = -f^sym " 6L, (37) 
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and this expression has been extensively used in the hterature to characterize the isospin 
dependence of the incompressibihty of asymmetric nuclear matter 12, Ss], 75, 76, 77, 78 1. 
Obviously, we have 



-f^sat,2 — -^"^asy TT^^ (38) 

and the coefficient i^asy thus could be a good approximation to -ft'sat,2 if the 3rd-order deriva- 
tive of the EOS of symmetric nuclear matter with respect to density, i.e., Jq, is negligible or 
the magnitude of the slope parameter L for the symmetry energy is very small. It should be 
noted that the higher-order derivatives of the EOS of symmetric nuclear matter with respect 
to density in Eq. (jl]), such as Jq, do not contribute to -R'sat,2- In the following, we will check 
how the Jo term affects the value of -ft'sat,2- 

As shown in Appendix „ the expressions Eq. (!23|) for the saturation density, Eq. (127|) 
for the binding energy and Eq. flM|) for the isobaric incompressibihty coefficient are exact 
up to 4th-order in 5. It is thus interesting to see that with precision up to 4th-order in (5, 
we only need to know 11 coefficients -Eo(po), Kq, Jo, h, -Ssym(po), L, K^y^, J^y^, -Esym,4(po), 

-^sym,4; -^'^sym,4 

among the 14 coefficients Eq{pq), Kq, Jq, Jq, Esymipo), L, Ksymi Jsjmi Isymi 
Esym,4.{po), -Z^sym,4, -f^sym,4, Jsym,4 and /sym,4 which are defined at the normal nuclear density 
Pq. The higher-order coefficient /sym in Eq. f|TT]) for the symmetry energy as well as Jsym,4 
and /sym,4 in Eq. fll6l) for the 4th-order symmetry energy do not affect the saturation density, 
the binding energy and the isobaric incompressibihty coefficient with precision up to 4th- 
order in 6. Furthermore, the 6 coefficients Eo{po), Kq, Jq, Esym{po), L, and i^sym determine 
completely the saturation density, the binding energy and the isobaric incompressibihty 
coefficient with precision up to 2nd-order in 6. 



III. THEORETICAL MODELS 



In this section, we introduce the theoretical models used in the present work and also 



give some important expressions for comp 
range Gogny effective interaction (MDI) 



eteness. These models include the modified finite- 



43( 1 . the Hartree-Fock approach based on Skyrme 



interactions, and a phenomenological modified Skyrme-like (MSL) model. A very useful fea- 
ture of these models is that analytical expressions for many interesting physical quantities 
in asymmetric nuclear matter at zero temperature can be obtained, and this makes it phys- 
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ically transparent and very convenient to check the higher-order effects on the properties of 
asymmetric nuclear matter. 



A. Isospin- and momentum-dependent MDI interaction 



The isospin- and momentum- dependent MDI interaction is based on the finite-range 



Gogny effective interaction 



431]. In the MDI interaction, the potential energy density V{p, 6) 



of an asymmetric nuclear matter at total density p and isospin asymmetry 6 is expressed as 



follows 



43 



55| 



-A-uPnPp -^l f 2 , 2\ 



B p' 



,(7+1 



a 



1 P^o 



Po 



X 




/ frir,p)fr'{r,f^ 



(39) 



i + {p-p')yA^' 

In the above, r = 1/2 (—1/2) for neutrons (protons); cr = 4/3; frif^p) is the phase-space 
distribution function of neutrons or protons at coordinate r and momentum p. The param- 
eters Au{x)^ Ai{x), B^Cr^T^Cr-T and A are obtained by fitting the momentum-dependence 
of the single-nucleon potential U (p, 5, p, r) to that predicted by the Gogny Hartree-Fock 
and/or the Brueckner-Hartree-Fock calculations, the saturation properties of symmetric nu- 
clear matter, and the symmetry energy at the normal nuclear matter density po 43|. The 
parameters Au{x) and Ai{x) are given by 

2B 



IB 

Au{x) = -95.98 - X -, Ai{x) 



-120.57 



X- 



(40) 



in terms of the parameter x, which is introduced to mimic various i?sym(p) predicted by 



different microscopic and/or phenomenological many-body theories [55|. By adjusting the x 
parameter, the i?sym(p) is varied without changing any property of symmetric nuclear matter 
and the symmetry energy at the saturation density Esym{po), as the x-dependent Au{x) and 
Ai{x) are adjusted accordingly. Using the definition in Eq. ([2]), we have -E'sym(Po) = 30.5 
MeV at Po = 0.16 fm"^ while its value is 31.6 MeV within the parabolic approximation of 
Esjm{p) ~ E{p,6 = 1)—E{p,6 = 0). We note that the MDI interaction has been extensively 
used in the transport model for studyin g is ospin effects in intermediate energy heavy- ion 



collisions induced by neutron-rich nuclei 



14 



15 



44 



46 



55 



56, 



79 



80 



81 



8^, the 



13 



study of the thermal properties of asymmetric nuclear matter 



85 



86 



density and pressure in neutron star crust 



72 



871 ] , and the transition 



73| . In particular, the isospin diffusion data 



from NSCL/MSU have constrained the value of x to between and —1 for nuclear densities 



from about 0.3po to 1.2po 



55, 



56|. 



With /r(^, p) = ■p-0(p/(r) — p) for nuclear matter at zero temperature, the integrals in 



Eq. (!39|) can be evaluated analytically. In particular, we have 
J J ^ ^ 1 + [p-p'Y/K^ 

+ 4A[(p^(r) - p){t')) arctan 



42 



vAr)-pf{T') 



ip^ij) + p){t')) arctan 



A 
A 



X In 



+ ^[A^ + QK\p){t) + p){t')) - Kpj{r) - p]{T')f] 



{Pf{r)+Pf{T')) 



iPfir)-Pf{T')y + A\ 
The kinetic energy per nucleon of cold asymmetric nuclear matter is 



(41) 



47r 



ip)n + P%)^ 



(42) 



hmh'^p 

where Pfn{p) = ^(3vr^p„(p))^/'^ is the Fermi momentum of neutrons (protons). The total 
energy per nucleon of cold asymmetric nuclear matter can be obtained as 



P 



(43) 



By setting Pn = Pp = 2 ^^^Pfn = Pfp = Pf, where = ?i(37r^p/2)^/^ is the fermi momentum 
of symmetric nuclear matter at density p, we obtain following EOS for cold symmetric nuclear 
matter: 



14 



Svr 



■p) + ^{Ai + Au) + 



B 



X {v){^v) - A^) - 8Apj arctan ^ 



o- + 1 VPo 



+ -[A^ + 12A2p5]ln 



4p2 + A2 



Furthermore, from Eq. fH3|) the symmetry energy can be expressed as 

-C/symlPj 



2V55V5=o 

S^mh?p'^ 4po ' " ^ + 1 VPo. 



+ 



+ 



9poP V ^ V 



4p) - A In 



4p^ + A^ 
A^ 



o o o 4n2. + A^ 
4p/-P/(4P/ + A')ln ^ 



9pQp \ J 

and the 4th-order nuclear symmetry energy can be written as 
1 fd^E 

Esyra AP) " 77 



A2 



4! V d5^ 



5=0 



Cu /47r 



7AV/ In 



4p2 + A2 4(7A^p^ + 42A2p6 + 40pS 



A2 



3%p V^^ 



A^ 



(4p2 + A2)^ 
A2 A2 



'(7AVJ + 16p)) In + 



(44) 



(45) 



(46) 



B. Skyrme-Hartree-Fock approach 



In the standard Skyrme Hartree-Fock model [6j, l88|, l89|, l90|, l91|, l92|, l93|, [94, l95|, l96| , the 
interaction is taken to have a zero-range, density- and momentum-dependent form with the 
interaction parameters determined from fitting the binding energies and charge radii of a 
large number of nuclei in the periodic table. In this approach, the EOS of asymmetric 
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nuclear matter can be expressed as [9l|, [92, l95|, [96 1 



10m V 2 



+ -top[2(xo + 2) - (2x0 + l)F2] 

o 

+ It3p-+l[2(x3 + 2) -(2x3+1)^2] 



+ 40 



3 /37r2\'/=^ 5 



where 



X {[ti(xi+2)+t2(x2 + 2)]F5/3 

+ i[t2(2a:2 + l)-ti(2a;i + l)]F8/3}, 



(47) 



The EOS of symmetric nuclear matter can thus be written as 



3h^ 




2/3 


lOm 


K 2 J 








> 2/3 


+ -Qs 




80 







2 

p3 



3 



lo 

with 0s = 3ti + (5 + 4x2)^2- Furthermore, the symmetry energy can be obtained as 

1 /d^E' 



(48) 



-E'sym(p) 



2 V ,=0 



6m \ 2 



2 1 . 

p:^ - -to(2Xo + l)p 



1 fsn^y/' , 

-2a[—) 

- 1^3(2X3 + l)p'^+\ 



(49) 



with 9sym = 3tiXi — ^2(4 + 5x2). Similarly, the 4th-order nuclear symmetry energy can be 
written as 



EsymAP) - 77 



1 /d^E\ 


4! V -95^ ; 






162m 






^37r' 


' 648 I 


^ 2 



<5=0 
2\ 2/3 



2 
p3 



2\ 2/3 



0sym,4P" 
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with 6sym,4 = - xi) + 3^2(1 + 3^2)- 



C. A phenomenological modified Skyrme-like (MSL) model 



Following the energy density functional obtained from the Hartree-Fock approach with 
the zero-range, density- and momentum-dependent form with the Skyrme interaction, the 
binding energy per nucleon of a cold asymmetric nuclear matter at total density p and isospin 
asymmetry 6 in the modified Skyrme-like (MSL) model is parameterized as 



Emsl{p, S) 



.5/3 
rn 



2m. 



-.pT 



a p 

+ 77- + - 

2 po 0- 



sym 



(50) 



where = | (Stt^)^^^; a, (3 and a are parameters; and E^^^^{p) represents the local density- 
dependent part of the symmetry energy. The effective neutron and proton masses m* and 



m* are assumed to have forms similar to the standard SHF results 



91 



92 



95 



96|, i.e.. 



2m* 
_^ 
2m* 



2m 
2m 



p (Ccff + /^cff^) 
P (Cefj - Doff^) 



(51) 
(52) 



where Ceg and Defr are constants. This implies that the single-nucleon potential depends 
quadratically on the nucleon momentum as in the standard SHF approach. In terms of the 
isoscalar effective mass m* and the isovector effective mass m* given by 

2m 

_^ 

2m 



2m* 



2m* 



+ pCcs 

+ p (Ceff - Deff) 



(53) 

(54) 



the nucleon effective mass can be written as 97 1 



^2 



2p„ 



2m* 



1 



2p, 
P 



2m* 



q = n,p. 



(55) 



„g p 2m* 

We note that the isovector effective mass m* corresponds to the proton (neutron) effective 
mass in pure neutron (proton) matter. Also, we can easily obtain the following relation 



2m* 2m!: 



26 



2m*. 2m* 



(56) 
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Experimentally, the isoscalar effective mass m* and the isovector effective mass m* at the 
normal nuclear density po hay e been constra i ned to be m*Q ^ 0.8m and m*Q ^ 0.7m, 



respectively 



42 



91 



92 



98 



99, 



IQO, 



101 



102 



103l |. With these constrained values for the 



isoscalar and isovector effective masses, Eq. ( l56l) gives a larger neutron effective mass than 
the proton effective mass in a neutron-rich matter, which is consistent with experimental 
data on the isospin dependence of the nucleon optical potential and also recent microscopic 
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33, 



45|. 



and phenomenological many-body theory predictions 

The EOS of symmetric nuclear matter in the MSL model is then given by 

2/3 / N 5/3 



Eo{p) 



-^kin 



PoJ 



C 



PoJ 



a p 
2 po (7 + 1 



(57) 



2/3 



2/3 

Po 



where the first term represents the kinetic energy contribution with E'^^^^ — 
and the second term is due to the nucleon effective mass with the coefficient C being a 
constant determined by the isoscalar effective mass m* q as 



C 



m 



m 



m 



(58) 



s,0 



The parameters a, (3 and a in the MSL model are determined by the binding energy 
per nucleon Eq{pq) and the incompressibility Kq of cold symmetric nuclear matter at the 
saturation density po? and they can be expressed as 

2 



a 



X 



2i?kV 



K-^kin 

- IOC 



3Eo(po) - 2C) 



X 



Ko + 9i?o(po)-i?L 

g EM - hj) 

Ko - 9Eo{po 



E^ 

•-^kin 



^-^kin 



4C 



16C 



a 



Ko + 9Eo{po)-E^,^-4C 
Kq + 2E°, - IOC 
3E,V-9Eo(po)-6C- 



(59) 

(60) 
(61) 
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In particular, for Eo{po) = —16 MeV, m*Q = 0.8m and po = 0.16 fm ^, we have 

C = 5.53 (MeV) (62) 

a = -47.90 - 39.37—^ — (MeV) (63) 

Kq — 188.21 

_ /fo- 11.05 

177.16 ' ^ ^ 

where the units of Kq is MeV. 

The symmetry energy in the MSL model can be expressed as 

/ \ 2/3 / \ 5/3 

EsyM = i?sym(P0) (j-J + ^sym(p), (66) 

2 / Q 2 \ 2/3 

where the first term is the kinetic energy contribution with E^^^{pq) = ^ ( 2 Po) ^^"^ 
the second term is due to the contribution of the nucleon effective mass with the coefficient 
D determined by both m* q and m* q as 

(67) 



9 V <,o <,o 



With m*Q = 0.8m and m*^ = 0.7m at po = 0.16 fm ^, we have D = —3.51 MeV. In Eq. fl66l) . 



similarly to the momentum- independent MID model 5l|, the local density-dependent part 



-^svm(p) is parameterized as 



sym 

'Ysym 



^srm(p) = (1 - yXm(Po)^ + yEl;i,ipo) ^ ) (68) 

Po VPo/ 

with the constant EI°'^{pq) determined by 

^s°ym(P0) = ^sym(Po) - E^^ipo) - D. (69) 

Obviously, we have ^'"^(po) = 21.2 MeV following Esym(Po) = 30 MeV and E^j^{po) = 12.3 
MeV at Po = 0.16 fm~^. The default value for the 7sym parameter is taken to be 4/3 in 
the MSL model following the -E'sym(p) in the MDI interaction, namely, Eq. 0451) (we will 
see how the '~fsym parameter affects the symmetry energy in the following). In particular, 
similarly to the x parameter in the MDI interaction, the dimensionless parameter y in the 
MSL model is introduced to mimic various i?sym(p) predicted by different microscopic and/or 
phenomenological many-body theories for a fixed 7sym parameter. As we will show later, for 
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7sym = 4/3, adjusting the y value can nicely reproduce the Esjm{p) in the MDI interaction 
with X = —1, 0, and 1. Therefore, the symmetry energy density functional constructed in 
the MSL model is very flexible and can mimic very different density behaviors by varying 
only one parameter. 

In the MSL model, similarly to the SHF approach, the 4th-order and higher-order nu- 
clear symmetry energies only include contributions from the kinetic energy and the nucleon 
effective mass while the local density-dependent part of higher-order symmetry energies are 
neglected in Eq. (150|) . In particular, the 4th-order nuclear symmetry energy in the MSL 
model can be shown to be 



where the first term is the kinetic energy contribution while the second term is due to the 
contribution from the nucleon effective mass. 



The MSL model is thus an extension of the momentum- independent MID model [5l|] by 
including the effects of the nucleon effective mass. It provides a simple phenomenological 
parametrization of the EOS of asymmetric nuclear matter and is thus a convenient and 
transparent way to investigate the possible correlations among higher-order and lower-order 
characteristic parameters of asymmetric nuclear matter. In the MSL model, we have totally 
8 free parameters, i.e., C, D, a, /3, a, Esym{po), y, and 7sym which can be determined 
by empirical information on the EOS of symmetric nuclear matter, the nucleon effective 
mass and the density dependence of symmetry energy. In particular, the parameters C and 
D (or equivalently Ccs and -Dcfr) are determined by the isoscalar effective masses m* and 
the isovector effective mass m* at the normal nuclear density po, i.e., m*Q and r«^o- The 
parameters a, (3 and a are determined by Eo{pq), Kq and po while the parameters Esjm{po), 
y, and 7sym are introduced to mimic the density dependence of different symmetry energies 
predicted by microscopic and/or phenomenological many-body theories. As a default in 
the MSL model, we use m*g = 0.8m, m* Q = 0.7m, po = 0.16 fm~^, -E'o(po) = -16 MeV, 
Kq = 240 MeV, -Esym(Po) = 30 MeV, and 7sym = 4/3 and vary the parameter y to describe 
different symmetry energies. 
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IV. RESULTS AND DISCUSSIONS 



A. Characteristic parameters at normal nuclear density and EOS of asymmetric 
nuclear matter 



As shown in Section [TTl the expressions for the saturation density Eq. ([2SD, the binding 
energy Eq. ( 1271) and the isobaric incompressibihty coefficient Eq. (IMl) are exact up to 
4th-order in 6, and these expressions involve 11 characteristic parameters defined at the 
normal nuclear density po, i-e., -Eo(po), Kq, Jq, Iq, E^ymipo), L, K^ym, Jsym, E^ymAPo), 
Lsym,4, and Ksym,A- We summarize the values of po, -Eo(po), Ko, Jo, h, Kasy, K^„t,2 and K^.^t,2 

in Tablesdlini andEIlwhile the values of Esym(po), L, K^ym, Jsym, -Esym,4(Po), Lsym,i, KsymA, 

m*gQ/m and m*g/m in Tabels HV] W\ and IVII for the MDI interaction with x = 1,0,-1 
and the popular 63 standard Skyrme interactions with their saturation density and the 
symmetry energy satisfying 0.140 fm~^ < po < 0.165 fm~^ and 25 MeV< -E'sym(Po) < 37 
MeV, respectively. For the 63 standard Skyrme interactions, the values in the Tables are 
sorted in the order of increasing values of L. It should be stressed here that the parameters 
of all Skyrme interactions are chosen to fit the binding energies and charge radii of a large 
number of nuclei in the periodic table. Detailed values of the parameters for these 63 Skyrme 

50, 89, 9^, 91, 92, 93, 95, 96 



interactions can be found in Refs. 



104 



105 



106 



107 



108 



110 



111 



112 



1091 . 



1131 ]. The selected ranges of po and -EsymlPo) are essentially consistent with 



their empirical values inferred from experimental data. We note that here no constraints are 
imposed on Kq and L for selecting the Skyrme interactions as we will systematically explore 
the correlations of other physical quantities with Kq or L. 

Since the 11 characteristic parameters are defined at the normal nuclear density po, it is of 
interest to study the extent these characteristic parameters can provide information on the 
properties of asymmetric nuclear at sub-saturation and supra-saturation density regions. As 
an example, we show in Fig. [T]the energy per nucleon of symmetry nuclear matter from the 
MDI interaction as a function of its density. Also shown in Fig. [T] are the results obtained 
by using Eq. (jlj) including terms up to x^, X^, and x'^, respectively. It is seen that Eq. 
(jl]) with terms up to x"^, the parabolic approximation Eq. flTUl) which involves only 
two characteristic parameters, i.e., -E'o(po) and Kq, can approximate very well the EOS of 
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TABLE I: The saturation density po and the characteristic parameters Eq{pq) (MeV), Kq (MeV), 
Jo (MeV), Iq (MeV), K^y (MeV), Ksat,2 (MeV) and Ksat,2 (MeV) at saturation density for the 
MDI interaction with x = —1,0,1 and the SHF predictions with 63 standard Skyrme interactions. 
The smah differences from Table I of Ref. [73| for some Skyrme interactions are due to the use of 

n 

0.17 (0.33) as an approximation of 1/6 (1/3) for the a parameter in Ref [731 1. The results shown 
here are thus more accurate than those in Ref [731 1. 
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-447.6 


2160.8 


-352. 


.0 


-321.1 


-8.4 


MDI(O) 0. 


160 


-16.2 


212, 


,5 


-447.6 


2160.8 


-443. 


.1 


-316.3 


52.9 


TV TT~V T / -1 \ 
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,160 
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212, 


,5 


-447.6 
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.3 
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Z 


0. 


159 
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,3 
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.6 
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0. 
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,6 


-352.4 
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.6 
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E 


0. 
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,5 
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,1 
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41.1 




0. 
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,3 


-369.0 


1546.0 


-225. 


.1 
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43.3 


SVI 


0. 


,143 
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,6 
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-1107.4 


-427. 


.3 


-424.2 


-4.4 


z* 


0. 
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,9 


-369.2 


1544.4 


-305. 


.5 


-312.6 


-7.1 


SkSC4 


0. 


,161 


-15.9 


234. 


,7 


-380.8 


1549.2 
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.5 


-320.1 


-4.1 


SI 


0. 


,155 
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,4 
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.2 


-469.7 


-14.2 


BSk3 


0. 


,157 


-15.8 


234. 


,8 


-380.9 
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,6 
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BSkl 


0. 


,157 
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,3 


-385.6 


1588.7 


-325. 


.0 


-313 


-8.4 


SIII 


0. 


,145 


-15.9 
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,4 


101.4 


-903.0 


-453. 


,2 


-456.0 


-20.2 


BSk2 


0. 


,157 


-15.8 


233. 


,7 


-380.1 


1542.4 


-344. 


,8 
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0. 
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,2 


-385.4 


1587.3 
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.0 


-315.4 
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0. 
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,8 


-367.2 
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,3 
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.0 
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0. 


,157 


-15.8 
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,1 


-370.6 


1571.3 


-316. 


.3 


-289.0 


-16.8 


BSk7 


0. 


,157 


-15.8 


229. 


,3 


-370.9 


1572.8 


-317. 


.3 


-288.2 


-17.0 


SKP 


0. 


163 
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201. 


,0 


-435.6 


2127.8 
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,3 


-341.9 


7.4 
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0. 


,157 
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237. 


,2 
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,8 
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,159 
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,1 


-380.4 
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-435. 
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,160 
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Force po Eo{po) Kq Jq 
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TABLE II: Continued with Table U 



Force 




Po 


^o(po) 


Ko 


Jo 


lo 


K 


Ksat,2 




SKX 


0, 


.155 


-16. 


,1 


271, 


,1 


-297.4 


904.0 


-451, 


,2 


-414. 


,8 


2.3 


SKXce 


0, 


,156 


-15. 


,9 


268. 


,2 


-294.6 


892.9 


-439, 


,3 


-402. 


,5 


3.8 


BSkl5 


0, 


.159 


-16. 


,0 


241, 


,6 


-363.1 


1457.0 


-395, 


,9 


-345, 


,4 


-2.7 


BSkl6 


0, 


.159 


-16. 


,1 


241, 


,7 


-363.6 


1459.9 


-396, 


.6 


-344, 


,2 


-1.7 


BSklO 


0, 


.159 


-15. 


,9 


238, 


,8 


-370.3 


1479.6 


-418, 


.3 


-360, 


,6 


12.7 


SGII 


0, 


.158 


-15. 


,6 


214, 


,7 


-380.9 


1741.8 


-371, 


,7 


-304, 


,9 


17.8 


BSkl2 


0, 


.159 


-15. 


,9 


238, 


,1 


-369.1 


1474.9 


-419, 


,4 


-360, 


,5 


14.0 


BSkll 


0, 


.159 


-15. 


,9 


238, 


,1 


-369.2 


1475.3 


-420, 


.0 


-360, 


,5 


14.6 


SLylO 


0, 


.156 


-15. 


,9 


229, 


,7 


-358.3 


1559.6 


-374, 


,7 


-314, 


,2 


-24.6 


BSkl3 


0, 


.159 


-15. 


,9 


238, 


,1 


-369.2 


1475.2 


-420, 


.8 


-360, 


,6 


15.5 


BSk9 


0, 


.159 


-15. 


,9 


231, 


,4 


-374.9 


1591.5 


-384, 


,7 


-320, 


,0 


-3.1 


BSkl4 


0, 


.159 


-15. 


,9 


239, 


,3 


-358.7 


1434.8 


-415, 


.5 


-349, 


,7 


14.2 


SLy230a 


0, 


.160 


-16. 


,0 


229, 


,9 


-364.2 


1593.6 


-364, 


,1 


-293, 


,9 


-32.1 


SLy6 


0, 


.159 


-15. 


,9 


229, 


,9 


-360.2 


1568.8 


-383, 


,7 


-312, 


,9 


-13.2 


SLy8 


0, 


.160 


-16. 


,0 


229, 


,9 


-363.2 


1587.1 


-388, 


.4 


-316, 


,8 


-12.0 


SLy4 


0, 


.160 


-16. 


,0 


229, 


,9 


-363.1 


1586.9 


-392, 


.1 


-320, 


,5 


-12.7 


SLyO 


0, 


.161 


-16, 


,0 


230, 


,2 


-365.2 


1598.7 


-389, 


.2 


-317, 


,2 


-11.6 


SLy3 


0, 


.160 


-16, 


,0 


229, 


,9 


-363.4 


1588.0 


-395, 


.4 


-323, 


,4 


-12.9 


SKM* 


0, 


,160 


-15. 


,8 


216, 


,6 


-386.1 


1768.9 


-430, 


.6 


-349. 


,0 


37.3 


SLy230b 


0, 


.160 


-16, 


,0 


229, 


,9 


-363.1 


1586.8 


-395, 


.5 
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,9 


-12.2 


SLy7 


0, 


.158 
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,9 
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,7 


-359.2 
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.7 


-327, 


,6 


-10.4 




n 
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1 Cr Q r r 
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- ( .u 
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0, 


.160 
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,8 


216, 


,6 


-386.1 
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.9 


-356, 


,9 


45.1 


SII 


0, 


.148 


-16, 


,0 


341, 


,4 


15.8 


-567.5 


-565, 


.9 


-568, 


,2 


23.2 


SLy5 


0, 


.161 


-16, 


,0 


229, 


,9 


-364.1 


1592.7 


-413, 


.7 


-334, 


,0 


-4.6 
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TABLE III: Continued with Table HI 



Force 




Po 


Eoipo) 


Ko 


Jo 


lo 


K 


Ksat,2 


Ksat,A 


SLy9 


0. 


,151 


-15.8 


229.8 


-350.4 


1511.4 


-413.7 


-329.2 


1.7 


SkI6 


0. 


,159 


-15.9 


248.2 


-326.6 


1251.8 


-402.2 


-324.3 


15.5 


SkI4 


0, 


,160 


-15.9 


248.0 


-331.2 


1280.0 


-402.9 


-322.2 


21.9 


SIV 


0, 


,151 


-16.0 


324.6 


-68.8 


-234.9 


-517.7 


-504.2 


64.8 


SGI 


0. 


,154 


-15.9 


261.8 


-297.9 


1005.0 


-435.2 


-362.5 


61.8 


SKO* 


0. 


,160 


-15.7 


222.1 


-390.2 


1706.7 


-495.6 


-373.2 


98.3 


SkMP 


0. 


,157 


-15.6 


231.0 


-338.4 


1425.2 


-468.9 


-366.8 


101.5 


SKa 


0, 


,155 


-16.0 


263.2 


-300.1 


1014.4 


-526.2 


-441.1 


105.0 


SKO 


0. 


,160 


-15.8 


222.8 


-391.3 


1712.3 


-519.1 


-379.5 


143.9 




0. 


,158 


-15.6 


237.4 


-348.4 


1377.2 


-523.3 


-397.5 


184.5 


Gcr 


0. 


,158 


-15.6 


237.2 


-348.7 


1379.5 


-550.1 


-411.9 


230.6 


SKT4 


0. 


,159 


-16.0 


235.5 


-383.0 


1562.5 


-589.3 


-436.2 


201.2 


SV 


0. 


,155 


-16.0 


305.7 


-175.8 


183.5 


-552.4 


-497.1 


176.2 


SkI3 


0. 


,158 


-16.0 


258.2 


-303.9 


1088.3 


-530.1 


-411.8 


152.1 


SkI2 


0. 


,158 


-15.8 


240.9 


-339.7 


1351.4 


-555.3 


-408.2 


226.9 


SkI5 


0. 


,156 


-15.8 


255.8 


-302.0 


1083.7 


-616.4 


-463.7 


347.2 



symmetric nuclear matter from about 0.5po to l.Spo- Including higher-order terms of 
and with the characteristic parameters Jq and Jq in Eq. (jlj) improves significantly the 
approximation to the EOS at low densities and that up to about 2po- To describe reasonably 
the EOS of symmetric nuclear matter above 2po, one needs to include higher-order terms in 
X ■ 

Fig. [2] displays the density dependence of symmetry energy using the MDI interaction 
with X = 1, 0, and —1 together with corresponding results obtained by using Eq. ( ITT|1 
up to X? X^i X^) s-iid respectively. It is seen that the importance of the contributions 
from higher-order terms in x to the density dependence of symmetry energy depends on the 
stiffness of symmetry energy. For a supra-soft symmetry energy (x = 1), terms up to x^ ^ire 
needed to describe reasonably the symmetry energy from sub-saturation densities to about 
2po- The situation is similar for the case of modestly soft symmetry energy (x = 0). For the 
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TABLE IV: The characteristic parameters -Esym(yOo) (MeV), L (MeV), Ksym (MeV), Jsym (MeV), 
-E-sym^C/Oo) (MeV), i^sym,4 (MeV), -fCsym,4 (MeV), 171* and m*Q/m at saturation density for the 
MDI interaction with j; = — 1, 0, 1 and the SHF predictions with 63 standard Skyrme interactions. 



Force 


_hj Qiim ( 0C\ ) 
byTil \rV } 


L 


K 

syvn 


sym 


Esym,4{Po) Lsym,A 






m* n/ui 


MDI(l) 


30.5 


14.7 


-264.0 


660.0 


U.oz 


U.oo 


-4.82 


0.67 


0.54 


MDI(O) 


30.5 


60.2 


-81.7 


295.3 


U.oz 


U.oo 


-4.82 


0.67 


0.54 


MDI(-l) 


30.5 


105.8 


100.7 


-69.3 


U.OZ 


U.Oo 


-4.82 


0.67 


0.54 


z 


26.8 


-49.7 


-657.9 


495.2 


n 

U. ( 


z.Oo 


1.47 


0.84 


0.73 




26.4 


-36.9 


-457.8 


880.0 


U.oo 


O.U4 


2.38 


0.84 


0.70 


E 


27.7 


-31.3 


-570.7 


448.6 


n en 


Z.OO 


1.67 


0.87 


0.74 




26.7 


-29.4 


-401.4 


883.1 


U.yi 


Q 1 /I 

0.i4 


2.58 


0.78 


0.66 


SVI 


26.9 


-7.3 


-471.3 


146.0 


U.o ( 


z.Uo 


0.76 


0.95 


0.81 


z* 


28.8 


-4.5 


-332.6 


725.1 


n no 


O.Z4 


2.79 


0.77 


0.65 


SkSC4 


28.8 


-2.2 


-329.5 


708.3 


U.4D 


U.yi 


-1.85 


1.00 


1.00 


SI 


29.2 


1.2 


-461.8 


141.4 


U. ( U 


O 1 

z.io 


0.73 


0.91 


0.80 


BSk3 


27.9 


6.8 


-306.9 


550.3 


n VI 
U. ( i 


z.iy 


0.76 


1.12 


0.89 


BSkl 


27.8 


7.2 


-281.8 


606.4 


n /I Q 


n vn 
U. ( y 


-2.04 


1.05 


1.05 


BSk2 


28.0 


8.0 


-297.0 


557.9 


n VI 
U. ( i 


Z.io 


0.74 


1.04 


0.86 


MSk7 


27.9 


9.4 


-274.6 


592.1 


U.4o 


n vo 
u. ( y 


-2.04 


1.05 


1.05 


SIII 


28.1 


9.9 


-393.7 


130.4 


U.oo 


O CO 

z.oy 


2.34 


0.76 


0.66 


BSk4 


28.0 


12.5 


-265.9 


558.4 


0.61 


1.70 


-0.22 


0.92 


0.85 


BSk8 


28.0 


14.9 


-220.9 


624.9 


0.43 


0.78 


-2.09 


0.80 


0.87 


BSk6 


28.0 


16.8 


-215.2 


603.5 


0.45 


0.89 


-1.85 


0.80 


0.86 


BSk7 


28.0 


18.0 


-209.4 


598.2 


0.42 


0.77 


-2.08 


0.80 


0.87 


SKP 


30.0 


19.6 


-266.8 


508.6 


0.94 


3.33 


2.96 


1.00 


0.74 


BSk5 


28.7 


21.4 


-240.3 


499.9 


0.64 


1.83 


0.04 


0.92 


0.84 


SKXm 


31.2 


32.1 


-242.8 


428.7 


0.88 


3.02 


2.40 


0.97 


0.75 


RATP 


29.2 


32.4 


-191.2 


440.6 


1.06 


3.94 


4.21 


0.67 


0.56 
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TABLE V: Continued with Table HV] 



Force 


Esym 


(Po) 


L 




Jsym 




Ksym,i 






SKX 


31 


.1 


33.2 


-252.1 


379.7 


0.89 


3.10 


2.61 


0.99 


0.75 


SKXce 


30, 


.1 


33.5 


-238.4 


356.9 


0.89 


3.12 


2.65 


1.01 


0.75 


BSkl5 


30 


.0 


33.6 


-194.4 


466.5 


0.64 


1.84 


0.03 


0.80 


0.77 


BSkl6 


30 


.0 


34.9 


-187.4 


461.9 


0.60 


1.66 


-0.31 


0.80 


0.78 


BSklO 


30 


.0 


37.2 


-194.9 


397.0 


0.69 


2.08 


0.50 


0.92 


0.81 


SGII 


26 


.8 


37.6 


-145.9 


330.4 


0.87 


3.01 


2.39 


0.79 


0.67 


BSkl2 


30 


.0 


38.0 


-191.4 


392.5 


0.68 


2.02 


0.40 


0.92 


0.82 


BSkll 


30 


.0 


38.4 


-189.8 


390.1 


0.67 


2.01 


0.37 


0.92 


0.82 


SLylO 


32 


.0 


38.7 


-142.2 


591.2 


0.37 


0.50 


-2.60 


0.68 


0.80 


BSkl3 


30 


.0 


38.8 


-187.9 


386.6 


0.67 


2.00 


0.37 


0.92 


0.82 


BSk9 


30 


.0 


39.9 


-145.3 


475.8 


0.37 


0.47 


-2.71 


0.80 


0.91 


BSkl4 


30 


.0 


43.9 


-152.0 


388.3 


0.60 


1.65 


-0.33 


0.80 


0.78 


SLy230a 


32 


.0 


44.3 


-98.2 


602.9 


0.06 


-1.07 


-5.80 


0.70 


1.00 


SLy6 


31 


.2 


45.2 


-112.5 


511.3 


0.38 


0.55 


-2.54 


0.69 


0.80 


SLy8 


31 


.4 


45.3 


-116.5 


511.4 


0.40 


0.64 


-2.39 


0.70 


0.80 


SLy4 


31, 


.8 


45.4 


-119.9 


521.0 


0.40 


0.63 


-2.40 


0.69 


0.80 


SLyO 


31, 


.5 


45.4 


-116.8 


510.6 


0.40 


0.64 


-2.39 


0.70 


0.80 


SLy3 


32 


.1 


45.5 


-122.1 


526.2 


0.40 


0.63 


-2.41 


0.70 


0.80 


SKM* 


30, 


.0 


45.8 


-155.9 


330.5 


0.94 


3.32 


2.97 


0.79 


0.65 


SLy230b 


32 


.0 


46.0 


-119.7 


521.5 


0.40 


0.61 


-2.43 


0.69 


0.80 


SLy7 


32 


.4 


48.1 


-114.3 


516.6 


0.38 


0.55 


-2.53 


0.69 


0.80 


SLy2 


32 


.3 


48.8 


-113.5 


502.9 


0.40 


0.63 


-2.39 


0.70 


0.80 


SLyl 


32 


.5 


48.8 


-115.7 


508.5 


0.40 


0.63 


-2.40 


0.70 


0.80 


SKM 


30 


.7 


49.3 


-148.8 


323.3 


0.91 


3.19 


2.71 


0.79 


0.66 


SII 


34, 


.2 


50.0 


-265.7 


104.7 


1.10 


4.21 


4.94 


0.58 


0.50 


SLy5 


32 


.7 


50.3 


-111.9 


499.2 


0.40 


0.63 


-2.40 


0.70 


0.80 
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TABLE VI: Continued with Table |V] 



Force 


EsymiPo) 


L 


K 

^ ^ sym 


J sym 


Esym,4:{Po) Lsym,A 


Ksym,A 






SLy9 


32.1 


55.4 


-81.3 


461.8 


0.33 


0.32 


-2.88 


0.67 


0.80 


SkI6 


29.9 


59.2 


-46.8 


378.1 


0.28 


0.04 


-3.55 


0.64 


0.80 


SkI4 


29.5 


60.4 


-40.6 


351.2 


0.30 


0.15 


-3.36 


0.65 


0.80 


SIV 


31.2 


63.5 


-136.7 


79.5 


1.37 


5.51 


7.50 


0.47 


0.41 


SGI 


28.3 


63.9 


-52.0 


194.5 


0.86 


3.00 


2.36 


0.61 


0.57 


SKO* 


32.1 


69.7 


-77.5 


221.4 


0.55 


1.38 


-0.90 


0.90 


0.87 


SkMP 


29.7 


69.8 


-50.3 


159.7 


0.93 


3.31 


3.00 


0.65 


0.58 


SKa 


32.9 


74.6 


-78.5 


174.5 


1.13 


4.33 


5.06 


0.61 


0.51 


SKO 


32.0 


79.6 


-42.3 


130.0 


0.59 


1.57 


-0.53 


0.89 


0.85 


Ro- 


30.6 


85.7 


-9.1 


22.2 


0.85 


2.88 


2.13 


0.78 


0.68 


G(j 


31.4 


94.0 


14.0 


-26.7 


0.85 


2.87 


2.12 


0.78 


0.68 


SKT4 


35.5 


94.1 


-24.5 


97.8 


0.45 


0.91 


-1.83 


1.00 


1.00 


SV 


32.8 


96.1 


24.2 


48.0 


1.70 


7.18 


10.77 


0.38 


0.33 


SkI3 


34.8 


100.5 


73.0 


211.5 


0.12 


-0.74 


-5.10 


0.58 


0.82 


SkI2 


33.4 


104.3 


70.7 


51.6 


0.37 


0.48 


-2.66 


0.68 


0.80 


SkI5 


36.6 


129.3 


159.6 


11.7 


0.12 


-0.72 


-5.04 


0.58 


0.80 



stiffer symmetry energy (a; = — 1), including terms up to already give a good description 
of the symmetry energy from sub- saturation densities to about 3po- 

In Fig. [31 we show the density dependence of the 4th-order symmetry energy -Esym,4(p) 
using the MDI interaction together with corresponding results obtained by using Eq. ( |T6l) up 
to X; X^; X") and x"^, respectively. We note here that the -Esym,4(p) does not depend on the x 
parameter in the MDI interaction as shown in Eq. P6l) . Firstly, one can see from Fig. [3]that 
the magnitude of -Esym^lp) is very small compared to that of Esymip)- As shown in Table HVl 
the value of -E'sym,4(p) at the normal nuclear density po is about 0.62 MeV which is consistent 
with the predictions from the SHF approach using different Skyrme interactions as shown 
in Tables HVt IVl and I VII where one can see that only 5 Skyrme interactions, i.e., RATP, SII, 
SIV, SKa, and SV, among the 63 Skyrme interactions have -E'sym.4(po) larger than 1 MeV 
(but still less than 2 MeV). The value of -Esym,4(Po) = 0.62 MeV is further consistent with 
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FIG. 1: (Color online) Energy per nucleon as a function of density for symmetric nuclear matter 
in the MDI interaction. Also included are results obtained by using Eq. ([3]) up to x^i X'^i 
respectively. 

the value 0.57 MeV predicted by the MSL model using Eq. ( 170|) . These results thus confirm 
the empirical parabolic law that the higher-order (including 4th-order) contributions of 6 in 
the EOS of asymmetric nuclear matter are usually very small and negligible as mentioned 
previously. Furthermore, similarly to the cases shown in Figs. [1] and [21 including terms up 
to ill Eq. (fT6|) can approximate very well the exact 4-th-order symmetry energy for the 
density region between about 0.5po and 1.5po while including higher-order terms of and 
X'^ improves significantly the results at low densities and up to about 2po- 

The above results thus indicate that generally one needs higher-order terms in x (higher 
than 4th-order) to describe the EOS of asymmetric nuclear matter at high density region 
(above 2po). We note that above conclusions obtained with the MDI interaction are also valid 
for the SHE approach and the MSL model. These features imply that it is very difficult to 
obtain correct information on the high density behaviors of the EOS for asymmetric nuclear 
matter based on the characteristic parameters obtained at the normal nuclear density po- At 
this point, it should be mentioned that the transport model analysis of heavy-ion collisions 
at intermediate and high energies as well as the astrophysical observations, especially on 
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FIG. 2: (Color online) Density dependence of symmetry energy using the MDI interaction with 
X = 1 (left) (middle), and —1 (right) together with corresponding results obtained by using Eq. 
pT]) up to X: X^: and respectively. 



the properties of compact stars, provide unique too^ 
asymmetric nuclear matter at high densities j^, 12 



s to extract information on the EOS of 



15 



26 



114|. 



B. Isospin dependence of the saturation properties of asymmetric nuclear matter 

In the following, we show the results on the saturation properties of asymmetric nuclear 
matter, i.e., the saturation density as well as the binding energy and incompressibility at 
saturation density. Especially, we investigate their isospin dependence and study if the 
higher-order terms in the isospin asymmetry 6 (5^ term) are important for the description 
of the saturation properties of asymmetric nuclear matter. 



1. The saturation density 

The saturation density is a basic quantity of asymmetric nuclear matter. To see the 
symmetry energy dependence of the saturation density of asymmetric nuclear matter, we 
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FIG. 3: (Color online) Density dependence of the 4th-order symmetry energy -E'sym,4(p) using 
the MDI interaction together with results obtained by using Eq. ()16p up to X: X^i X'^j ^-iid 
respectively. 

use here the MDI interaction with a; = 1, 0, and —1. The density dependence of the 
symmetry energy from this interaction is shown in Fig. |H Also included in Fig. [ 1] are 
the results from the widely used APR (Akmal-Pandharipande-Ravenhall) prediction 115 1 
and the phenomenological MSL model prediction with y = —15, —7.5, and 0.75 (We will 
discuss the MSL results later) . It is seen that the APR prediction for the symmetry energy 
resembles very well that from the MDI interaction with x = up to about 3.5po- 

Using the MDI interaction with x = 1, 0, and —1, we have calculated the density and 
isospin asymmetry dependence of the binding energy per nucleon of asymmetric nuclear 
matter, and the results are shown in Fig. [51 Further indicated in Fig. [5] are corresponding 
saturation points in the E-p plane. One can see that different symmetry energies lead to 
different EOS of pure neutron matter (5 = 1) as expected. In particular, the EOS of pure 
neutron matter for x = —1 is bounded at low densities. In addition, different symmetry 
energies lead to rather different behaviors for the saturation points in the E-p plane. 

In order to see more clearly the isospin dependence of the saturation density, we show 
in Fig. [6] the saturation density Psat(<^) as a function of 6"^ in the MDI interaction with 



30 




FIG. 4: (Color online) Density dependence of the symmetry energy from the MDI interaction 
with X = 1, 0, a nd —1. The results from the widely used APR (Akmal-Pandharipande-Ravenhall) 



prediction 



115( 1 and the phenomenological MSL model prediction with y = —15, —7.5, and 0.75 



are also included for comparison. 

X = 1, 0, and —1, respectively. The exact saturation density Psat(^) is obtained from Eqs. 
( !2T|) and (143|) . Corresponding results from Eq. (l23l) including terms up to 5^ and up to 
6^, respectively, are also included for comparison. The results indicate that the saturation 
density generally decreases with isospin asymmetry and more neutron-rich nuclear matter 
has lower saturation density. In addition, for the stiffer symmetry energy (x = —1), the 
nuclear matter can be bounded even for pure neutron matter (The corresponding saturation 
density is about 0.3po)- The inset in Fig. [6] displays corresponding results at smaller isospin 
asymmetries with 6"^ < 0.1 which is relevant to the properties of finite nuclei. In the 
small isospin asymmetry region (5^ < 0.1), the saturation density Psat(5) displays a linear 
dependence on 6"^ and Eq. fl2^ including terms up to 6"^ thus approximates very well the 
exact Psat(^)- Furthermore, how Psat(^) decreases with increasing 6 depends on the stiffness 
of symmetry energy with the softer symmetry energy having a weaker dependence, and this 
feature can be easily understood from Eq. flM|) which indicates that the slope of Psat(<^) 
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FIG. 5: (Color online) The density and isospin asymmetry dependence of the binding energy 
per nucleon for asymmetric nuclear matter in the MDI interaction with x = 1, 0, and —1. The 
saturation points at different isospin asymmetries are also indicated. 

with respect to 5^ is proportional to — ^po- In addition, at larger isospin asymmetries with 
5^ > 0.3, including higher-order 6 terms up to 5^ in Eq. ( l23l) still deviates significantly 
from the exact /Osat(^) and higher-order terms of S are thus necessary (except for the case of 
X = where Eq. fl23l) with terms up to 6"^ gives a good approximation to the exact Psat(<^) 
in the whole 6 region where the asymmetric matter can still be bounded). These features 
imply that higher-order terms in 6 may be important for the determination of the saturation 
density of nuclear matter at very neutron-rich nuclear environment, such as inside a neutron 
star. 

2. Binding energy at saturation density 

The isospin dependence of the binding energy per nucleon of asymmetric nuclear matter 
at saturation density Esati^) is shown in Fig. [7] as a function of 5^ for the MDI interaction 
with X = 1,0, and —1. The exact -Esat(<^) is obtained from Eqs. fl2T|) and fH3|) . Corresponding 
results from Eq. (p7j) including terms up to 6"^ and up to 6^, respectively, are also included 
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FIG. 6: (Color online) The saturation density psat as a function of 5^ in the MDI interaction with 
X = 1, 0, and —1. Corresponding results from Eq. ()23p including terms up to 5^ and up to 5^, 
respectively, are also included for comparison. The inset in Fig. [6] displays corresponding results 
at smaller isospin asymmetries with 6'^ < 0.1. 

for comparison. It is seen that i?sat(^) generally increases with increasing isospin asymmetry. 
The results at smaller isospin asymmetries with 5^ < 0.1 are shown in the inset of Fig. [7] 
and it is seen that i?sat(5) displays a linear dependence on 6^ and therefore can be very 
well approximated by Eq. fl27j) including terms up to 6"^. We note that the rate at which 
-^'satl^) increases with 6"^ at small 6"^ is determined uniquely by -E'sym(Po) a-s shown in Eq. 
( 127|) . Also, it is seen that including higher-order 6 terms up to 5^ in Eq. ( 1271) gives a good 
approximation to the exact Esa,t{S) in the whole 6 region (except for the case of x = —1 
where Eq. (1271) including terms up to 6^ still deviates significantly from the exact -E'sat(^) at 
larger isospin asymmetries with 6"^ > 0.3 and higher-order terms of 6 are thus needed). 

3. Incompressibility at saturation density 

Shown in Fig. [8] is the incompressibility at saturation density -K'sat(^) as a function of 
5^ for the MDI interaction with x = 1, 0, and —1. The exact i^'sat(5) is obtained from 
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FIG. 7: (Color online) Same as Fig. [6] but for the binding energy at saturation density i^sat- 

Eqs. (1211) . fl52]) . and fH51) . and corresponding results from Eq. flMl) including terms up to 
5^ and up to 5^ are also included for comparison. It is seen that Jisat(^) generally decreases 
with increasing isospin asymmetry and more neutron-rich nuclear matter has smaller in- 
compressibility. This feature is consistent with earlier calculations based on microscopic 



many-body approaches 1161 ]. The softening of the incompressibility of asymmetric nuclear 



matter with increasing isospin asymmetry may have imp ortant implications in understand- 



ing the mechanism for supernovae explosions |116l . Ill7l |. Corresponding results at smaller 
isospin asymmetries with 5^ < 0.1 are given in the inset of Fig. [HI and it shows that Eq. 
f lM|) including terms up to 5^ approximates very well the exact A'sat(<^) as Ksa.t{S) displays 
a good linear correlation with 5^. As to the decreasing rate of i^sat(^) with increasing 5^ 
at small 6"^, it is determined by the parameter -ft'sat,2 as shown in Eq. flM|) which depends 
on the characteristic parameters Jo, Kq, L, and Ksym- In addition, including higher-order 
6 terms up to 5^ in Eq. ( l34l) is seen to give a good approximation to the exact Ksat{S) in 
the whole 6 region (except for the case of x = —1 where Eq. (1341) including terms up to 
5^ deviates significantly from the exact Ksati^) at larger isospin asymmetries with 5^ > 0.3 
and higher-order terms of 6 are thus important). 

The above results indicate that the saturation properties of asymmetric nuclear matter, 
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FIG. 8: (Color online) Same as Fig. [6] but for the incompressibility at saturation density i^^sat- 

i.e., the saturation density as well as the binding energy and incompressibility at saturation 
density, exhibit a good linear correlation with 5^ at smaller isospin asymmetries with 6^ < 0.1 
which is relevant to the properties of finite heavy nuclei. On the other hand, depending on the 
stiffness of nuclear symmetry energy, higher-order terms in 6 {6^ and higher-order terms) may 
become important for describing reasonably the saturation properties of asymmetric nuclear 
matter at larger isospin asymmetries with 6^ > 0.3. The importance of higher-order isospin 
asymmetry terms for the stiffer symmetry energy has also been observed in previous studies 



on the transition density in neutron stars 



72 



731 ]. In addition, the saturation density and 



the incompressibility at saturation density generally decrease with the magnitude of isospin 
asymmetry while the binding energy at saturation density shows an opposite behavior. 
Again, we note that above conclusions obtained from the MDI interaction are also valid for 



the SHF approach and the MSL mode 



recent study based on the RMF model |ll8 |. 



Our results are further consistent with the very 
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C. Constraining the Ksat,2 parameter from the phenomenological MSL model 

1. General information on the -ftrsat,2 parameter 
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FIG. 9: (Color online) The absolute values of i^sat,2 and Ki^y and the value of -fCsat,4 as functions 
of L for the MDI interaction with x = 1, 0, — 1 and the 63 standard Skyrme interactions considered 
in the present work. 



As shown in the above, the i<'sat,2 parameter essentially characterizes the isospin de- 
pendence of the incompressibility of asymmetric nuclear matter at saturation density and 
the higher-order parameters (e.g., -R'sat,4) are only important for extremely neutron-rich (or 
proton-rich) nuclear matter with stiffer nuclear symmetry energies. Actually, it can be seen 
from Tables [H [TTl and IIIII that the magnitude (absolute values) of -ft'sat,2 is generally much 
larger than that of -ft'sat,4 for the MDI interaction with x = 1,0,-1 and the 63 standard 
Skyrme interactions considered in the present work. Shown in Fig. [9] are the absolute values 
of -ft'sat,2 and i^asy as well as the value of -R'sat,4 as functions of L for the MDI interaction with 
X = 1, 0, — 1 and the 63 standard Skyrme interactions considered in the present work. It is 
seen that these values of -R'sat,2 can be nicely expressed as —400 ± 120 MeV. For the magni- 
tude of the i^asy parameter, it is generally larger than that of the -ft'sat,2 parameter, especially 
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for the stiffer symmetry energies (larger L values), which indicates that the higher-order Jq 
parameter is important as discussed later. Furthermore, the absolute values of -R'sat,2 are 
clearly much larger than that of -ft'sat,4 except that at very large L values the absolute values 
of -R'sat,4 may become larger and comparable with that of -R'sat,2- This feature is consistent 
with the results shown in Fig. [S] where the higher-order terms are seen to be only important 
for the stiffer symmetry energies. 

It is generally believed that information on - K'sa.t.2 can be extracted experimentally by 



measuring the GMR in neutron- rich nuclei 119l |. Usually, one can define a finite nucleus 



incompressibility Ka{N, Z) for a nucleus with neutrons and Z protons [A = N + Z) hy 
the energy of GMR -Egmr, i-e., 



_ l h^KA{N,Z) 

- V m (r2) ' ^^^^ 

where m is the nucleon mass and (r^) is the mean square mass radius of the nucleus in the 
ground state. Similarly to the semi-e mpir ical mass formula, the finite nucleus incompress- 



ibility Ka{N, Z) can be expanded as 119 1 



2 



+i^Coui^ + ■ ■ •, (72) 

Neglecting the i^curv term, the K^s term and other higher-order terms in Eq. (1721) . one can 
express the finite nucleus incompressibility Ka{N, Z) as 

2 



(N — Z\' 



^2 

+Kcoui^^, (73) 

where Kq, i^surf, K^, and i^coui represent the volume, surface, symmetry, and Coulomb 
terms, respectively. The parameter is usually thought to be equivalent to the -ft'sat,2 
parameter. However, we would like to stress here that the -ft'sat,2 parameter is a theoretically 
well-defined physical property of asymmetric nuclear matter as shown previously while the 
value of the pa rame ter may depend on the details of the truncation scheme in Eq. (!72l) . 
As shown in Ref. jl20| . K^. may be related to the isospin-dependent part of the surface 
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properties of finite nuclei, especially the surface symmetry energy. Therefore, cautions are 
needed to interpret the parameter as the Ksa.t,2 parameter (we will go back to this point 
later) . 

Earlier attempts based on the above method have given widely different values for the 
Kj. parameter. F or ex ample, a value of Kj. = —320 ± 180 MeV with a large uncertainty was 
obtained in Ref. 12 1| from a systematic study of the GMR in the isotopic chains of Sn and 



Sm. In this analysis, the value of Kq was found to be 300 ± 25 MeV, which is somewhat 
larger than the commonly accepted value of 240 ± 20 MeV. In a later study, an even less 
stringent constraint of —566 ± 1350 < Kr < 139 ± 1617 MeV was extracted from the GMR 
of finite nuclei, depending on the mass region of nuclei and the number of parameters used 



122]. More recently, a much more 



in parameterizing the incompressibility of finite nuclei 
stringent constraint of AV = —550 ± 100 MeV has been obtained in Ref. 
measurements of the isotopic dependence of the GMR in even-A Sn isotopes. 



22 



231] from 



2. Correlation between Jq and Kq 



As shown in Eq. fl35l) . the -R'sat.2 parameter is completely determined by the 4 char- 
acteristic parameters Kq, Jq, L, and -K'sym at the normal nuclear density. It thus would 
be interesting to estimate the possible value of -ft'sat,2 from knowledge on Kq, Jq, L, and 
i^sym- For the incompressibility of symmetric nuclear matter at its saturation density po, 
the transport model analyses on experimental data from subt hreshold kaon production in 



heavy-ion collisions favor a soft equation of state 



26|, 



123 



1241 ]. More recently, the value of 



Kn has 
y, 21 



Deen more stringently determined to be 240 ± 20 MeV from the nuclear GMR data 



22 



23 



24]. 



While the Kq parameter has been relatively well determined, the Jq parameter is poorly 
known, and there is actually no experimental information on the Jq parameter. In the MSL 
model, one can easily calculate from Eq. fl57j) the Jq parameter as 



27p; 



1 



\p=po 



Ku. - ^Eq 



-\Kl 



2C 

+12EI^Eq + 6C(3E0i, - 25Eo)]. 



18^0 - 2E°, - SC)Kq 



(74) 



38 



> 

03 



200 


-200 

-400 
1 



- 



■1 - 



1 1 1 1 1 1 1 1 
_ ■ SHF 


1 1 1 1 1 


1 1 1 1 1 


^ MDI 

yj\ IVIUI 






" MSL (m *=m) 

^ s,0 ' 






~ . ivisL (m *=0.9m) 






" MSL (m^J=0.8m) 






~ MSL (m^J=0.7m) 












— — \ — \ — \ — \ — — \ — \— 

- 


H — \ — — \ — h- 
1 1 1 1 1 


H — \ — — \ — \ — 

- 

1 1 1 1 1 


200 250 


300 


350 



\ (MeV) 

FIG. 10: (Color online) The third derivative parameter Jq and its ratio Jq/Kq to the incom- 
pressibility Kq of symmetry nuclear matter at saturation density as functions of Kq from the MSL 
model, the MDI interaction and the SHF prediction with the 63 Skyrme interactions shown in 
Tables Ulnl andlffll 



In particular, for m* q = 0.8m, we have 

Jo = [Kl - 376.4A'o + 11214.645) (MeV), 

while we have 

Jo = W - 332.2/^0 - 4243.2) (MeV), 

for m* Q = m (and then C = 0) and 

1 



Jn 



51.2 



[Kl - 408^0 + 22252) (MeV) 



(75) 



(76) 



(77) 



for m*g = 0.7m and C = 9.47 MeV. In above equations, the unit of Kq is MeV. Therefore, 
the Jo parameter in the phenomenological MSL model is a quadratic function of the Kq 
parameter. 
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Shown in Fig. [10] are Jo and Jq/Kq as functions of Ko in the MSL model for different 
values of m*Q = m, 0.9m, 0.8m and 0.7m. Also included in Fig. [10] are corresponding 
results from the MDI interaction and the SHF prediction with the 63 Skyrme interactions 
shown in Tables [T] [ITl and IIIII It is seen that the correlation between Jq and Kq is similar 
among these three different models. Also, all three models show an approximately linear 
correlation between Jq/Kq and Kq. This linear correlation can be easily understood from 
Eqs. f[75|) . f[76|) and f[77|) . For example, on the r.h.s of Eq. f[75]) . the last term is small 
compared with the first two terms, and thus one has Jq ~ {Kq — 376 AKq), and then 
Jq/Kq ^ ^-y {Kq — 376.4) with the units of Kq in MeV. The linear dependence of Jq/Kq 
on Kq becomes better for larger values of m* q which usually leads to smaller values for the 
last two terms on the r.h.s of Eq. ( 1741) as shown in Eq. ( 1761) . We note that the correlation 
between Jq a nd Kq obtained in the present work is also consistent with the early finding by 



Pearson 



125|. In addition, it is seen from Fig. [TO] that the m*Q only has visible effects on 



the correlation between Jq and Jq/Kq for smaller Kq values while it has almost no effects for 
larger Kq values. While there do not exist any empirical constraints on the Jq parameter, 
we assume in the present study the correlation between Jq and Kq from the MSL model is 
valid and then determine Jq/Kq from the experimental constraint on Kq. 



3. Correlation between L and K, 



sym 



The parameters L and Ksym are determined by the density dependence of the symmetry 
energy around saturation density. In recent years, significant progress has been made both 
experimentally and theoretically in extracting the information on the behaviors of nuclear 
symmetry energy at sub-saturation density. Using the isospin- and momentum-dependent 
IBUU04 transport model with in-medium NN cross sections, the isospin diffusion data were 
found to be consistent with the symmetry energy from the MDI interaction with x between 
and —1, which can be parameterized by a density-dependent symmetry energy of -Esym(p) 
31.6(p/po)'^ with 7 = 0.69 — 1.05 at subnormal density (p < po) 



50 



55 



56|, leading to the 



extraction of 61 MeV < L < 111 MeV and -82 MeV < K^ym < 101 MeV 



50 



55, 



56| 



Using the Skyrme interactions consistent with the EOS obtained from the MDI interaction 



with X between and —1, the neutron-skin thickness of heavy nuclei calcu 



Hartree-Fock approach is consistent with available experimental data 50 



ated within the 



126| and also 
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FIG. 11: (Color online) Correlation between Xgym and L fr om the ]VISL model with ^sym — 

4/3 

and 5/3, the MDI interaction and the SHF prediction with the 63 Skyrme interactions shown in 
Tables IIV[ IVl and IVII The results from two simple one-parameter symmetry energies in Eqs. (jSip 
and (j82p are also shown for comparison. 



that from a relativistic mean-field model based on an accurately calibrated parameter set 
that repr oduc es the GMR in ^''Zr and ^°^Pb as well as the isovector giant dipole resonance 



of ^"«Pb 

-E'sym(p) 



127|. The extracted symmetry energy further agrees with the symmetry energy 
= 31.6(p/po)°'^^ recently obtained from the isoscaling analyses of isotope ratios 



65 MeV and Kgym ~ 



in intermediate energy heavy ion collisions [57[, which gives L 
—61 MeV. Furthermore, the above limited range of E^^^i^p) at subsaturation density is 
essentially consistent with the symmetry energy Esym{p) = 12.5(p/po)^^'^ + ^^-^ip/ PoV with 
7 = 0.4 — 1.05, extracted very recently from the ImQMD (Improved QMD) model analyses 



of both the isospin diffusion data and the double neutron/proton ratio 



The symmetry 



energy Esymip) = 12.5(p/po)^/^ + 17.6(p/po)'^ with 7 = 0.4 — 1.05 leads to the constraints of 
46 MeV < L < 80 MeV and -63 MeV < ifsym < -17 MeV. 

It should be noted that all above constraints on L and -R'sym are based on some unique en- 
ergy density functionals and thus special correlation between L and i^sym has been implicitly 
assumed. It is thus of interest to study if there exists a universal correlation between L and 
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Ksym in different models. For the MDI interaction, it is implied from Eq. (H5|) that the L 
and Ksyra both change linearly with the parameter x. Therefore they are linearly correlated 
by varying the parameter x that changes the density dependence of the symmetry energy. 
For the MSL model, we can obtain from Eq. fl^B]) following expressions 

L = 2E':^M + + 3 [i?sym(po) - i?^;m(Po) - D] 

+ 3y(7sym - 1) [E.yM ' i^^syZ(Po) ' D] (78) 
i^sym = 9?/7sym(7sym - 1) [^sym(po) - E^^iPo) " D] 

+ 10D-2E^;i{po). (79) 

Therefore, for the MSL model, Eq. (ITS]) and Eq. (I7U1) show that the L and -ft'sym both 
change linearly with the parameter y, and thus they are also linearly correlated by varying 
the parameter y to change the density dependence of the symmetry energy. In particular, 
we have 

i^sym = 37symi^ + ^sym(P0)(37sym-2) 

+2D{5 - 37sym) - 97sym^sym(Po)- (80) 

Also, L and -ft'sym are expected to be correlated within the SHE energy density functional. 
Shown in Fig. [TT] are the correlation between -ft'sym and L from the MSL model with 
= 4/3 and 5/3 (We used here the default values of ^sym(Po) = 12-3 MeV, ^sym(po) = 30 
MeV, and D = -3.51 MeV from ?Ti* g — 0.8?Ti and Tn*Q — 0.7?77,), the MDI interaction and 
the SHF prediction with the 63 Skyrme interactions shown in Tables IIV] IVl and IVI[ It 
is seen that the -R'sym parameter indeed displays approximately a linear correlation with 
the L parameter for the SHF predictions with the 63 Skyrme interactions and this linear 
correlation is nicely reproduced by the MDI interaction and the MSL model with '^sym = 4/3. 
For the MSL model, one can see from Eq. flHOl) that the 7sym parameter controls the shape 
(slope) of the linear correlation between L and -ft'sym- Furthermore, Fig. [11] shows that 
results from some Skyrme interactions deviate from the linear correlation obtained by the 
MDI interaction and the MSL model with 7sym = 4/3. In order to consider the uncertainty 
in the shape (slope) for the correlation between L and -ft'sym; we also include the result 
from the MSL model with '-/sym = 5/3. The correlation between -ft'sym and L from the SHF 
predictions with the 63 Skyrme interactions is consistent nicely with that from the MSL 
model with 7syjn = 4/3 and 5/3. Furthermore, we find that the isospin-dependent nucleon 
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effective mass has very small effects on the correlation between J'^sym and L. This can be 
easily understood from Eq. flHOl) since the value of 2D (5 — 37sym) is only about 7 MeV for 
7sym = 4/3 (it is zero for 7sym = 5/3). The linear correlation between -K'sym and L implies 



that one can extract -ft'sym from L. 



As mentioned above, an one-parameter parametrization for the symmetry energy is some- 
times used for simplicity, i.e.. 



E,y„(p) = Esym(po) (^J^ (81) 



or 

2/3 



i?sym(p) = £^Sm(P0) 

+ [i?sym(P0) - i?s';m(P0)] {£j . (82) 

For the parametrization of Eq. flSTl) . the L and K^ym can be obtained as 

L = 37Dsym(Po) (83) 
Ksjra = 97(7 - l)Esym(po, ) 



and furthermore we have 



^ -E'sym(po) 



which indicates that the -ft'sym parameter is quadratically correlated with L. The result 
obtained by varying the parameter 7 is shown in Fig. [H] by the dotted line. Similarly, for 
the parametrization of Eq. (!82|) . the L and -ft'sym can be expressed as 

L = 2E^^Jpo) + 37 [E^yM - i?syZ(Po)] (84) 

i^sym = 97(7 - 1) [EsyM ' E^^M] 

-2i?gS,(po), (85) 

and the i^sym parameter is again quadratically correlated with L, and this is illustrated by 
the dash-dot-dotted line in Fig. [TT] obtained by varying the value of 7. 

It is very interesting to see from Fig. [Ulthat for larger L values {L > 45 MeV which is 
consistent with the constraint from heavy- ion collision data shown later), all above symmetry 
energy functionals from different models and parameterizations give consistent predictions 
for the Ksy^-L correlation. This nice feature implies that using these different models and 
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parameterizations for the symmetry energy will not influence signiflcantly the determina- 
tion of the -ft'sat,2 parameter. On the other hand, the Ksym-L correlation from the two 
one-parameter parameterizations on the symmetry energies in Eqs. fl8ip and (152]) deviate 
signiflcantly from the MDI, MSL and SHF predictions for small L values. Actually, the two 
forms of one-parameter parametrization for the symmetry energy in Eqs. (1811) and (I82p may 
be too simple to describe a softer symmetry energy. As shown in Ref. 55|, Eq. (IHTl) cannot 
describe correctly the density dependence of the symmetry energy from the MDI interaction 
with X = 1 (the Gogny interaction). On the other hand, as shown in Fig. HI the MSL 
model can give a nice description of the density dependence of the symmetry energy from 
the very soft to the very stiff. Although there is no direct empirical information on the -ft'sym 
parameter and some uncertainty on the K^y^-L correlation still exist, we assume here that 
the correlation between i^sym and L from the MSL model with 7sym = 4/3 and 5/3 is valid 
and then use the experimental constraint on L to extract the value of -ft'sym- 



4- Phenomenological MSL model constraint on the i^sat,2 parameter 



From above analyses, we can now extract information on the -R'sat,2 parameter from the 
experimental constraints on the Kq parameter and the L parameter within the phenomeno- 



logical MSL model. As pointed out previously, the value of Kn has 



determined to be 240 ± 20 MeV from the nuclear GMR 



16 



21 



22 



Deen relatively well 



23 



2^. The slope 



parameter L has been found to correlate linearly with the neutron-skin thickness of heavy 



nuclei and thus can in princ i ple be de t ermi ned 



of such nuclei 



50 



64 



126 



128 



129 



130 



131 



rom measured thickness of the neutron skin 



132| . Unfortunately, because of the large 



uncertainties in the experimental measurements, this has not yet been possible so far. The 
proposed experiment of parity- violating electron scattering from ^°^Pb, i.e.. Parity Radius 
Experiment (PREx) at the Jefferson Laboratory is expected to give an indepe ndent and 
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134| and thus 



accurate measurement of its neutron skin thickness (within 0.05 fm) 
to impose a stringent constraint on the slope parameter L in future. On the other hand, 
as mentioned previously, heavy-ion collisions, especially those induced by neutron-rich nu- 
clei, provide a unique tool to explore the density dependence of the symmetry energy and 
thus the L parameter. Actually, as discussed previously, the L parameter has been already 
limited signiflcantly by heavy-ion collisions data. 
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FIG. 12: (Color online) ii'sat,2 as a function of L from the MSL model with 7sym = 4/3 (a) and 5/3 
(b) and m* Q = 0.8m and m* q = 0.7m for different values of Kq and -Esym(/Oo)- The shaded region 
indicates constraints within the MSL model with Kq = 240 ± 20 MeV, -Bsym(po) = 30 ± 5 MeV, 
and 46 MeV < < 111 MeV limited by the heavy-ion collision data. The results from the SHF 
approach with 63 Skyrme interactions are also included for comparison. In addition, the constraint 
of Kr = -550 ± 100 MeV obtained in Ref. 
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from measurements of the isotopic dependence 



of the GMR in even-A Sn isotopes is also indicated. 



From Eqs. fl5B|) and fISUl) . we have in the MSL model 



K 



sat, 2 



K 



sym 



+ 6-37; 



sym 



9Tsym-^sym (Po) • 



(86) 



Results on -ft'sat,2 as a function of L are shown in Fig. [12] for 7sym = 4/3 (panel (a)) and 5/3 
(panel (b)) and with Kq = 220, 240, and 260 MeV. In Eq. the value of Jq is obtained 

from Eq. (1711) using the value of Kq. For other quantities in the MSL model, the default 
(po) = 12.3 MeV, ^sym(po) = 30 MeV, and D = -3.51 MeV from m*Q = 0.8m 



values Egy^ 



and m*Q = 0.7m have been used in the calculations. Since the value of 2D{5 — 37, 



sym > 



IS 



only about 7 MeV for 7. 



sym 



4/3 and is zero for 7, 



sym 



5/3, as mentioned before, the 
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contribution of the D parameter to A'sat,2 is thus smalL Also, the correlation between -ft'sym 
and L depends on -Esymlpo) and D as shown in Eq. flHOj) . To take into consideration of the 
uncertainty in the value of -E'sym(Po)) we also include in Fig. [T2]the results with Kq = 220 
MeV and ^sym(po) = 25 MeV as well as Kq = 260 MeV and ^sym(Po) = 35 MeV, which 
represent, respectively, the upper and lower bounds for a fixed value of L. The shaded region 
in Fig. [T2] further shows the constrained L values from heavy- ion collision data, namely, 46 
MeV < L < 111 MeV. The lower limit of L = 46 MeV is obtained from the lower bound 
in the ImQMD analyses of the isospin diffusion data and the double neutron/proton ratio 
jsol l while the upper limit of L = 111 MeV corresponds to the upper bound of L from the 



IBUU04 transport model analysis of the isospin diffusion data 50|, l5J, |55|, |56| . The constraint 



46 MeV < L < 111 MeV is consistent with the analyses of the pygmy dipole resonances 



661], the giant dipole resonance (GDR) of ^ospb analyzed with Skyrme interactions 65|, the 



very precise Thomas- Fermi model fit to the binding energies of 1654 nuclei 63|], and the 



recent neutron-skin analysis [78|]. These empirically extracted values of L represent the best 
and most stringent phenomenological constraints available so far on the nuclear symmetry 
energy at sub-saturation densities. 

It is seen from Fig. [T2]that the -ft'sat,2 parameter decreases with increasing L for '~fsym = 4/3 
while the opposite behavior is observed for 7sym = 5/3. This feature can be easily understood 
from Eq. fl86l) . For %y-^ = 4/3, Eq. (186|) is reduced to 

Ks.t,2 = -(^ + 2)L - 12E,,M + 17.6 (MeV), (87) 

while for 7sym = 5/3, it is reduced to 

Ks.t,2 = -(^ + l)L - 15E,y„(po) + 36.9 (MeV). (88) 

For Ko = 240 ± 20 MeV, Fig. [H (or Eq. ([71])) shows that the value of Jq/^o varies from 
about —1.9 to —1.2 when m*Q changes from m to 0.7m. Eq. (IHTj) (Eq. (l88l) ) thus shows 
that Ksat,2 decreases (increases) with increasing L for jsym = 4/3 (5/3). 

An interesting feature observed from Fig. [T2] is that the K^^x.,! parameter depends signif- 
icantly on the symmetry energy at the normal nuclear density £'sym(po)- This can be seen 
more clearly from Eqs. ( IHTll and (l88ll which show that changing -E'sym(po) by 5 MeV leads 
to a variation of 60 — 75 MeV for Ksat,2- This feature indicates that an accurate determina- 
tion of -Esym(Po) is important for determining the value of Kss,t,2- In addition, we have used 
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Esymipo) = 30 MeV in Fig. [TT]and one can easily see that a variation of 5 MeV for i?sym(po) 
would lead to a shift of 60 — 75 MeV for the i^sym parameter for a fixed L. Therefore, the 
MSL model with 7sym = 4/3 and 5/3 together with -Esym(Po) = 30 ± 5 MeV may provide a 
nice estimate for both the K^y^-L correlation and its uncertainty. 

From the shaded region indicated in Fig. [T21 we find that for 7sym = 4/3, we have —437 
MeV < Ksat,2 < -292 MeV for L = 46 MeV while -487 MeV < K^^t,2 < -306 MeV for 
L = 111 MeV. For 7sym = 5/3, we have -477 MeV < K,^t,2 < -302 MeV for L = 46 MeV 
while -461 MeV < Ksat,2 < -251 MeV for L = 111 MeV. These results indicate that based 
on the MSL model with 4/3 < 7sym < 5/3, Kq = 240 ± 20 MeV, 25 MeV < ^sym(po) < 35 
MeV, and 46 MeV < L < 111 MeV, the A"sat,2 parameter can vary from -251 MeV to -487 
MeV. The results shown in Fig. [T2] is obtained from a Jq/Kq value that is evaluated with 
the default value m*Q = 0.8m. Similar analyses indicate that the Ksat,2 parameter varies 
from —261 MeV to —489 MeV if we use m*Q = 0.7m while it varies from —245 MeV to 
-485 MeV if TTi* Q — 0.9?Ti is used. These results indicate that the extracted value for /'i"sat,2 
is not sensitive to the variation of the nucleon effective mass. The MSL model analyses with 
4/3 < 7sym < 5/3, Kq = 240 ± 20 MeV, ^sym(po) = 30 ± 5 MeV, 46 MeV < L < 111 MeV, 
and m*Q = 0.8 ± 0.1m thus lead to an estimate of -R'sat,2 = —370 ± 120 MeV. 

Also included in Fig. [12] are the results from the SHF approach with the 63 Skyrme 
interactions. It is seen that among the 63 Skyrme interactions shown in Tables Ht IIT| and llllj 
there are only 19 Skyrme interactions with predicted -R'sat,2 consistent with the constraints 
obtained here (the shaded regions in Fig. [12] (a) or (b)). The 19 Skyrme interactions are 
SKa, SKT4, G,, Skl3, Skl2, R,, SKO, SKO*, SkMP, SGI, SKM, SLy5, SLyl, SLy9, SLy2, 
SLy7, SkI6, SLy230b, and SkI4. From Tables HI - IVI| one can see that all these 19 Skyrme 
interactions predict 0.151 fm'^ < po < 0.161 fm'^, -16.0 MeV < ^o(po) < -15.6 MeV, 217 
MeV < Ko < 263 MeV, 28 MeV < ^sym(po) < 35 MeV, and 46 MeV < L < 104 MeV, which 
are consistent with the empirical information. However, the value of the nucleon effective 
mass predicted by the 19 Skyrme interactions varies widely, i.e., from 0.58m to m for m*Q 
while from 0.51m to m for m* q. If we further impose the constraints of m*Q = 0.8 ± 0.1m 
and m*Q > m*Q, then only five Skyrme interactions, namely. Go-, Ro-, SKM, SKO and SKO* 
remain to be consistent with all known empirical constraints except that the SKM interaction 
gives a little smaller incompressibility Kq = 216.6 MeV than the empirical constrained values 
of Kq = 240 ± 20 MeV. These features imply that the constraint on the -ft'sat,2 parameter 
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does not significantly limit the nucleon effective mass as the former has previously been 
shown to be insensitive to the latter. 



5. Discussions on parameters i^sat,2) -f^asy, CL^d 

As shown in Eq. (l38l) . the K^sy parameter corresponds to the -R'sat,2 parameter when Jq is 
zero, i.e., the parabolic approximation to the EOS of symmetric nuclear matter (Eq. flTUl) ) 
is valid. In the MSL model, a vanishing Jg requires a Kq value of about 340 MeV, which 
is significantly larger than the empirical value of 240 ± 20 MeV. Since Jq/Kq ^ —1.5 for 
Kq = 240 MeV in the MSL model, we have -ft'sat,2 ~ -f^asy + or Kasy ~ -K'sat,2 — 1.5L. 
Therefore, the difference between the -ft'asy parameter and the -ft'sat,2 parameter depends on 
the slope parameter L of the symmetry energy and a larger L value (stiffer symmetry energy) 
would lead to a larger difference. It should be stressed, however, that the -ft'asy parameter 
is completely determined by the density dependence of symmetry energy regardless of the 
EOS of symmetric nuclear matter. 

Based on the IBUU04 transport model analysis of the isospin diffusion data, a value of 
-^asy = —500 ± 50 MeV has been extracted from the symmetry energy obtained by the 
MDI interaction with the x parameter between and —1. The extracted -ft'asy = —500 ± 
50 MeV is essentially consistent with the symmetry energy extracted from the ImQMD 
model analyses of both the isospin diffusion data and the double neutron/proton ratio jsol ]. 
which predict a -ft'asy value from about —500 MeV to —340 MeV, and the symmetry energy 



Esyiaip) = 31.6(p/po)°'^^ obtained from the isoscaling analyses |57|, which gives i^asy ~ —453 
MeV. Furthermore, the constraint i^asy = —500 ± 50 MeV is consistent with the very 
recent constraint of -ft'asy ~ — 500tioo MeV obtained from the study of neutron skin of 



finite nuclei 



7i 



]. In the MDI interaction, one finds Jq/Kq = —2.1 from Table [T] and we 
thus have —311 MeV < -ft'sat,2 < —316 MeV, corresponding to the prediction from the MDI 
interaction with the x parameter between and —1. Therefore, the -ft'sat,2 value for the 
MDI interaction is significantly larger than the -ft'asy value and is insensitive to the density 
dependence of symmetry energy. On the other hand, the Jq parameter in the MDI interaction 
is important for the determination of the -R'sat,2 parameter. The importance of the higher- 
order Jq parameter to A'"sat,2 can be further seen from Fig. [HI which shows that for the 
63 Skyrme interactions considered here, the magnitude of the i^asy parameter is generally 
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larger than that of the /'Csat,2 parameter, especially, for the stiffer symmetry energies (larger 
L values). These results indicate that the higher-order Jq contribution to Ksat,2 generally 
cannot be neglected. 



In Fig. [121 the constraint K^- = —550 ± 100 MeV obtained in Ref. 22|, |23|] from recent 
measurements of the isotopic dependence of the GMR in even-A Sn isotopes is also shown. 
While the estimate of /'^sat,2 = —370 ± 120 MeV obtained in the present work has small 
overlap with the constraint of Kr = —550 ±100 MeV, the latter still has significantly larger 
magnitude than the former. According to Fig. [121 there are only 6 Skyrme interactions, 
namely, SII, SIV, SV, SI, SkI5, SIII that predict J'S:sat,2 = -550 ± 100 MeV. As shown in 
Tables HI [TTl and llllt the values for Kq from SI, SII, SIII SIV, and SV are all larger than 305 
MeV which are obviously inconsistent with the empirical constraint of Kq = 240 ± 20 MeV. 
For SkI5, we have L = 129.3 MeV which is significantly larger than the empirical constraint 
of 46 MeV < L < 111 MeV. Therefore, none of the 63 Skyrme interactions considered in the 
present work is consistent simultaneously with -ft'sat,2 = —550 ±100 MeV and the empirical 
constraints of Kq = 240 ± 20 MeV and 46 MeV < L< 111 MeV. These features imply that 
the Kr parameter extracted from Eq. ( 1731) based on the GMR may not fully reflect the /'^sat,2 
parameter. As mentioned before, the parameter may depend on the detailed truncation 
scheme in Eq. ( !72|) . The constraint Kr = —550 ± 100 MeV obtained in Ref. 22|, |23| is based 



on Eq. f[73|) and thus neglects contributions from the K^urv term, Kgg and other higher-order 
terms in Eq. (|72l) . It is expected that including contributions from higher-orde r ter ms in 
Eq. (1731) may change the extracted value for the Kr parameter as found in Ref. |l22 |. 



V. SUMMARY AND CONCLUSIONS 



We have studied in the present paper higher-order effects on the properties of isospin 
asymmetric nuclear matter when the EOS and saturation properties of asymmetric nuclear 
matter are expanded in powers of the isospin asymmetry parameter 6 = {pn — Pp) / p and 
the dimensionless variable X = (p ^ Po)/3po that characterizes the deviations of the den- 
sity from the normal nuclear density po- Analytical expressions for the saturation density 
of asymmetric nuclear matter as well as its binding energy and incompressibility at satura- 
tion density have been derived and given exactly up to Ath-order in the isospin asymmetry 
^ = {Pn — Pp)/ P using 11 characteristic parameters defined at the normal nuclear density po 
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by the density derivatives of the binding energy per nucleon of symmetric nuclear matter, 
the symmetry energy Esym_{p) and the 4th-order symmetry energy i?sym,4(p), namely, -E'o(po), 
Kq, Jo, Iq, -Esym(Po), L, J-^sym, Jsym, -Esym,4 (po) , -^^sym,4, and KsymA' O^r method/recipc to 
derive the analytical expressions for the saturation properties of asymmetry nuclear matter 
is quite general, and in principle, the higher-order coefficients (higher than 5^) in the isospin 
asymmetry 6 can be easily obtained. Using the isospin- and momentum-dependent MDI 
interaction, the SHF approach with 63 popular Skyrme interactions, and the phenomeno- 
logical MSL model, we have systematically studied the higher-order effects on the properties 
of asymmetric nuclear matter. 

Firstly, our results indicate that including terms up to can give a good description 
of the EOS of symmetric nuclear matter -E'o(p), the density-dependent symmetry energy 
-E'sym(p), and the density-dependent 4th-order symmetry energy -E'sym,4(p) for densities less 
than about 2po while higher-order terms in x (higher than 4th-order and thus more charac- 
teristic parameters defined at po are necessary) are generally needed to describe reasonably 
the -E'o(p), -E'sym(p), and -E'sym,4(p) at higher density region (above 2po). These features imply 
that it is very difficult to predict the high density behaviors of the EOS of asymmetric nu- 
clear matter based on the characteristic parameters obtained at the normal nuclear density 
Po. Therefore, it is of great interest to use the transport model analysis of heavy- ion collisions 
at intermediate and high energies as well as the astrophysical observations of compact stars 
as tools to extract information on the EOS of asymmetric nuclear matter at high densities. 

Secondly, we have systematically studied the isospin dependence of the saturation prop- 
erties of asymmetric nuclear matter, i.e., the saturation density psat(5) as well as the binding 
energy Esax{5) and incompressibility i^'sat(^) at saturation density. In particular, we have 
compared their exact results with the analytical expressions expanded up to the 5^ term and 
the (5^ term, respectively. Our results indicate that the saturation properties of asymmetric 
nuclear matter, i.e., Psat(^), -E'sat(^) and Ksa.t{5), exhibit a good linear dependence on 5^ at 
smaller isospin asymmetries with 5^ < 0.1 which is relevant to the structure of finite heavy 
nuclei. This feature implies that the higher-order terms (the 5^ term and higher-order terms 
in 6) are not important for the saturation properties of asymmetric nuclear matter, at least 
for an asymmetric nuclear matter that is not extremely neutron-rich (proton-rich) . On the 
other hand, for asymmetric nuclear matter with extremely large isospin asymmetries with 
6"^ > 0.3 and depending on the stiffness of the nuclear symmetry energy, the higher-order 
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terms in S (5^ and higher-order terms) may become important for describing reasonably the 
saturation properties. It is further found that the saturation density and the incompressibil- 
ity at saturation density generally decrease with the magnitude of the isospin asymmetry 
while the binding energy at saturation density shows the opposite behavior. 

Finally, we have studied in detail the second-order isospin coefficient i^sat,2 of the incom- 
pressibility of an asymmetric nuclear matter at its saturation density. It is found that the 
magnitude of the higher-order -R'sat,4 parameter is generally small compared to that of the 
Ksa.t,2 parameter, so the latter essentially characterizes the isospin dependence of the incom- 
pressibility at saturation density. Furthermore, we have found that the -ft'sat,2 parameter is 
uniquely determined by L, -ft'sym and Jq/Kq. Since there is no experimental information on 
the Jo parameter and the -R'sym parameter, we have thus used the MSL model, which can 
reasonably describe the general properties of symmetric nuclear matter and the symmetry 
energy predicted by both the MDI model and the SHF approach, to estimate the value of 
-f^sat,2- Our results indicate that generally the higher- order Jq contribution to can- 
not he neglected, especially for larger L values. Interestingly, it is found that there exists a 
nicely linear correlation between Ksym (^nd L as well as between Jq/Kq one? Kq for the three 
different models used here, i.e., the MDI interaction, the MSL interaction, and the SHF 
approach with 63 Skyrme interactions. For the MSL model, the correlation between -R'sym 
and L is further found to depend significantly on the value of -Esym(Po) but not on variations 
of the nucleon effective mass. These correlations and features have allowed us to extract 
the values of the Jo parameter and the Ksym parameter from the empirical information on 
Ko, L and Esym{po)- In particular, using the empirical constraints of Kq = 240 ± 20 MeV, 
Esymipo) = 30 ± 5 MeV, A6 MeV < L < 111 MeV and m*o = 0.8 ± 0.1m m the MSL model 
leads to an estimate of Ksat,2 = —370 ± 120 MeV. 

While the estimated value of Ksat,2 = —370 ± 120 MeV in the present work has small 



overlap with the constraint of Kr = —550 ± 100 MeV obtained in Refs. [22, l23|] from recent 
measurements of the isotopic dependence of the GMR in even-A Sn isotopes, its magnitude 
is significantly smaller than that of the constrained Kr- Recently, there are several studies 



118 



135 



136| on extracting the value of the -R'sat,2 parameter based on the idea initiated 



by Blaizot and collaborators that the values of both Kq and -ft'sat,2 should be extracted from 
the same consistent theoretical model that successfully reproduces the experimental GMR 
energies of a variety of nuclei. These studies show that no single model (interaction) can 
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simultaneously describe correctly the recent measurements of the isotopic dependence of the 
GMR in even-A Sn isotopes and the GMR data of ^°Zr and ^osp]-, nuclei, and this makes 
it difficult to accura tely determine the value of -ft'sat,2 from these experimental data. Also, 



a very recent study 



137l | indicates that the effect due to the nuclear superfluidity ma y als o 



affect the extraction of the -R'sat,2 parameter from the nuclear GMR. As pointed out in 118| , 



these features suggest that the Kr = —550 ± 100 MeV obtained in Ref. [2^, l23| may suffer 
from the same ambiguities already encountered in earlier attempts to extract the Kq and 
Ksa.t,2 of infinite nuclear matter from finite-nuclei extrapolations. This problem remains an 
open challenge, and both experimental and theoretical insights are needed in future studies. 
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APPENDIX: DERIVATION OF ANALYTICAL EXPRESSIONS FOR THE SAT- 
URATION DENSITY AND INCOMPRESSIBILITY TO HIGHER-ORDER IN 
ISOSPIN ASYMMETRY 

For completeness, we present in this appendix a simple derivation of the analytical ex- 
pressions for the saturation density of asymmetric nuclear matter and its incompressibility 
at saturation density up to higher-order terms in the isospin asymmetry 6 = {pn — Pp) / p- 
In particular, analytical expressions for the saturation density and incompressibility at sat- 
uration density are given exactly up to 4th-order in 6. 

The binding energy per nucleon e(p, 6) of asymmetric nuclear matter can be expanded 
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up to 2nth-order in S as 

e(p, 5) = eo(p) + e2{p)5^ + e,{p)5'' + ■■■ + eMS'"". (A.l) 

The first term eo(p) represents the EOS of symmetry nuclear matter with 6 = 0. We further 
expand each ej(p) in Eq. (lA.ip up to nth-order in a dimensionless variable z = {p — po)/po, 
i.e. 

eo(p) = eo(po) + aoiz + ao2^^ H h flon^" 

e2(p) = e2(po) + flii^; + auz'^ H h ai„z" 

e^ip) = e^ipo) + (321-2 + 022-2^ H h a2„z" (A. 2) 

e2n(p) = e2„,(po) + ttniz + a„22^ H h a„n2;''. 

The coefficients aij represent the characteristic parameters of asymmetric nuclear matter 
defined at the normal nuclear density po- The saturation density psat is then determined by 
the following equation 



— 'p=P.at=0 (A.3) 



dp 

or equivalently 

||.=...=0 (A.4) 

with Zsat = (Psat - Po)/pO 

Substituting Eq. (lA.ip and Eq. (1A.2P into Eq. (1A.4P leads to the following equation 

aoi + 2ao22 + Boos^:^ H h nao„2;"'"^ + 

+ (aii + 2ai22 + 3ai32;^ H h nain2;""^)5^+ 

+ (a2i + 20222; + 3a23;2^ + ■ ■ ■ + na2nZ''-^)S^+ (A.5) 
+ ••■ + 

+ (a„i + 2a„22; + 3a„3z2 H h ?7.a„n2;"-"i)(5^" = 0. 

Obviously, we have Oqi = as the saturation density of symmetry nuclear matter is de- 



termined by |^|5=o,p=po = 0- order to obtain the solution to Eq. (lA.Sp . i.e., Zsa.t with 
precision of 0(5^'^"''^), we assume 



^sat = + ^4(5' + ■ ■ ■ + A2fc5'% (A.6) 

where (fn = 1, 2, fc) are coefficients to be determined from comparing the coefficients 
of 6"^^ terms (j = 1, 2, m) on the two sides of Eq. (lA.Sp after substituting Eq. flA.6p into 
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Eq. (lA.Sp . For example, for k = 1, i.e., ^sat = ^2^^, the coefficient A2 can be determined by 
the following equality 

2ao2^2 

+aii = 0. (A.7) 
Similarly, for k = 2, i.e., Zgat = ^26"^ + A^S^, the coefficient A4 can be determined by 

2ao2^4 + 3ao3^2 
+2012^2 

+021 = 0. (A.8) 

Furthermore, for /c = 3 and 4, the coefficients Aq and Ag can be determined, respectively, 
by 

2ao2 A + 6ao3^2 A + 4ao4^2 
+2012^4 + SaisAl 

+2a22^2 

+a3i = (A.9) 

and 

2002^8 + 8003(2^2^6 + Al) + UaoiAjA^ + 5ao5v4^ 
+2012^6 + 6013^2^4 + 4ai4A^ 
+2022^4 + 8023^2 
+2032^2 

+041=0. (A. 10) 

From above analyses, the coefficient of the 5^*" term, A2m, can be determined by all lower- 
order coefficients A2, A^^, A2m-2 ijn = 1, 2, 8, . . . , n). The A2m obtained above are complete 
and precise up to the order of 5^"*, and the higher-order characteristic parameters a^j do not 
contribute to A2m- In order to write down the general expressions for A2m, we define the 
following symbol 

Bn,m = A2 Al ■ ■ ■ A2n, (^-H) 
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where we have ii,i2, ■ ■ - in G {0, A^+} and they satisfy following conditions 



n. 



Then the coefficient of the S term satisfies the equation 

m— 1 m—j+1 
j=0 j=2 



from which we can obtain the coefficient as 

• m+l 

O02 



2m 



m—1 m—j+1 \ 

-\- ^ ^ ^ ^ "^OijiBi—i^rn—j ~l" '^ml j 
7 = 1 i=2 ^ 



(A.12) 



(A. 13) 



(A.14) 



Therefore, one can obtain the saturation density of asymmetry nuclear matter to any order 
of 5. 

Taking the precision to the order of 5^, we then have 



A2 
A, 



_ ail 

2ao2' 

ai2ail 



2a, 



3ao3af]^ 



a2l 



8aS, 2ao2 ' 



Converting the coefficients aj, into the conventional forms, i.e.. 



a02 - a03 - Y62' °04 - 24^' 



flu — -3 , ai2 



ai3 



6^21 — 3 ; 0,22 - 



18 ' 162 ' 

sym,4 



18 



we then obtain 



sym,4 



-"■0 







So the saturation density can be obtained as 

3L 



Psat — PO 



+ 



-"■0 



5^ + 0(5^ 



which is exactly Eq. ([22 



(A.15) 



(A. 16) 



(A.17) 



(A. 18) 
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The incompressibility coefficient Ksat of asymmetry nuclear matter at the saturation 
density is defined as 

Substituting Eq. 
obtain 

3-C/svm.4 3JqL^ 



Qp2 IP-Psaf 

T]) and Eq. (1A.2I) into Eq. (lA.lQjl and using Eq. 



(A.19) 

.181) . we can easily 



sat 




(A.20) 

which is exactly Eq. (15^ . We no te that the expression Ksat = Kq+ ^-ft'sym — 6L — 5^ 
was originally given in Ref. llQl ]. For terms higher than 5^, they can be straightforwardly 
obtained following above derivation but many more higher-order characteristic parameters 



aij would be needed. 
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